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A COMPARISON OF THE THURSTONE METHOD OF MULTI- 
PLE FACTORS WITH THE HOTELLING METHOD 
OF PRINCIPAL COMPONENTS 


C. H. McCLoy, ELEANOR METHENY, AND VIRGINIA KNOTT 
Iowa Child Welfare Research Station 


The factors found by the Hotelling method of principal compo- 
nents present the same necessity for rotation as those found by the 
Thurstone method of multiple factors. 

The Thurstone method of multiple factors is here found su- 
perior to the Hotelling method of principal components, and is sug- 
gested as the method of choice for all problems requiring the deter- 
mination of elemental components from a group of related variables, 
since it was found on every point of comparison that, when applied 
to a correlation matrix of unknown communalities, the Thurstone 
method was more accurate than the Hotelling method, and also re- 
quired much less time. 


Introduction 


In recent years two techniques have been presented for the pur- 
pose of extracting the elemental components from a group of related 
variables. In 1933 Thurstone* published A Simplified Multiple Fac- 
tor Method and an Outline of the Computations based upon what he 
called the “center of gravity” method. He subsequently revised this 
technique before presenting it in Vectors of Mind} in 1935. In 1933, 
independently of Thurstone, Hotellingt published “An Analysis of a 
Complex of Statistical Variables into Principal Components,” which 
was designed for the same purpose as the Thurstone multiple factor 
method. Thurstone early recognized the necessity of rotating his pri- 
mary factor loadings to give logically meaningful results, but Ho- 
telling has presented the results of his method as giving principal 
components without mentioning the need for rotation. However, as 
we shall demonstrate below, Hotelling’s factors should and may be 
rotated in exactly the same manner as Thurstone’s; so that we may 
state that the two techniques are but different means of arriving at 


*L. L. Thurstone, A Simplified Multiple Factor Method and an Outline of 
the Computations. (Autographic reproduction of typewritten copy) Chicago: 
The University of Chicago, 1933. 

+L. L. Thurstone, Vectors of Mind. Chicago: University of Chicago Press, 


t Harold Hotelling, “Analysis of a Complex of Statistical Variables into 
Principal Components.” Journal of Educational Psychology, XXIV: 417-441, 498- 
520. September and October, 1933. 


— 
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similar ends. It therefore becomes of interest to those using these 
methods to know their relative merits and shortcomings. 

Both methods present one difficulty, namely, that since the diago- 
nals of the correlation matrix are unknown, some approximation 
method must be devised to supply reasonably accurate diagonal val- 
ues. Theoretically, the values that should be inserted are those ob- 
tained by summing the squares of all of the “true” factor loadings, or 
what Thurstone has called the “communalities.” Since these are un- 
known, some substitute must be found. Thurstone inserts the largest 
value found in the row or column in each successive correlation or 
residual matrix. Hotelling uses the value 1.00 in all diagonals and 
retains the computed values on subsequent sheets. While in an analy- 
sis based on a large number of variables the error introduced by such 
approximations tends to be absorbed, in an analysis based on a small 
number of variables, this error may become quite large. Recognizing 
this as a source of error, the procedures suggested by the two authors 
were used in their respective methods in this study. 


Procedure 


To secure a valid criterion with which to compare the results of 
our analyses, points were located with reference to orthegonal axes 
X, Y, and Z as follows: 


Point Vector Location 
1 9 on X axis 
2 8 on X axis 
3 7 on Y axis 
4 8 on Y axis 
5 9 on Z axis 
6 6 on Z axis 
7 6 in plane YOZ equidistant from OY and OZ 
8 8 in plane XOY equidistant from OX and OY 
9 7 in plane XOZ equidistant from OX and OZ 
10 9 equidistant from planes XOZ, YOZ, and XOY 


giving the factorial matrix shown in Table I. 

Using the formula r = h, h. cos a, where h, and h, represent the 
vectors and a represents the angle between them, correlations were 
computed between the variables, giving the correlation matrix shown 
in Table II. 

The Thurstone center of gravity method was first applied and 
three factors were found. The absolute value of the residuals remain- 
ing after the third factor was extracted averaged .023. These unro- 
tated factor loadings are shown in Table III. 
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These factor loadings were rotated by successive approximations, 
according to Thurstone’s directions, and the rotated factor loadings 
are shown in Table IV. 

Using the same correlation matrix, Hotelling’s method was ap- 
plied, and again three factors were found, the average value of the 
third factor residuals being .096. These factor loadings are shown in 
Table V. 

The Hotelling factors were also rotated by the method of suc- 
cessive approximations, and these rotated values are shown in Table 
VI. 

The communalities found by the two methods were computea 
and compared with the original communalities. These values are 
shown in Table VII. 

Finally, for the purpose of comparing the mechanics of the two 
techniques independently of the errors introduced by the approxima- 
tions of the diagonal values, the postulated communalities were placed 
in the diagonals and the analyses repeated. In each analysis three 
factors were found, and in each the third factor residuals practically 
disappeared, those remaining being so small (not larger than .002) 
as to be accounted for by the successive rounding off of the residuals 
in computation. Both sets of factor loadings were then rotated by the 
method of successive approximations, and after rotation both were 
found to correspond almost exactly with the postulated factorial ma- 


trix. 
Discussion 

The correlations used are hypothetical and may be considered as 
exact and having no probable error. Since three factors were postu- 
lated, the computations should produce three factors, and the third 
factor residuals should disappear. 

It is noted that in neither method do the third factor residuals 
entirely disappear. This is to be expected, since in neither case, as 
was explained above, do the diagonal entries exactly correspond with 
the true communalities. Hence, even though the techniques are other- 
wise perfect, these necessary approximations would introduce a cer- 
tain amount of error which would be evidenced both in the factor 
loadings and in the residuals. The Thurstone residuals are smaller, 
the average of their absolute values being .023, as compared with 
.096 for the average of the Hotelling residuals. This would tend to 
indicate that Thurstone’s method more nearly approximates the true 
communalities. 

Turning to the unrotated factor loadings in Tables III and V, we 
note that neither set bears any relation to the original factor matrix. 
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This is in accord with Thurstone’s statement of the necessity for ro- 
tating the factors. 

After rotating both sets of factors, we note in Tables IV and VI 
that they now approach more or less closely the values found in the 
original factorial matrix. Since Hotelling’s values after rotation ap- 
proach the expected values of the original factorial matrix, whereas 
before rotation they bore no relation to those values, it would seem 
that to follow Hotelling’s directions and to attempt to interpret the 
unrotated factor values would be a fruitless proceeding. 

Noting the difference between the rotated factor values and the 
expected values, we find that the average error by Thurstone’s meth- 
od is .025, while that by Hotelling’s method is .054. 

One further point of comparison lies in the computed commu- 
nalities found by the two techniques, as presented in Table VII. It 
will be noted that the average of the differences between Hotelling’s 
communalities and the true values is .125, and that of Thurstone’s 
values is .056, which would tend to indicate that Thurstone’s method 
of approximating the communalities is more nearly correct than Ho-: 
telling’s. 

Thus it would seem that at every point of comparison the Thur- 
stone method shows itself superior in accuracy to the Hotelling 
method. 

One further point of superiority is also worthy of note. The com- 
putations by the Thurstone method may be completed very rapidly, 
while the iterative solution of Hotelling is not only tedious but ex- 
tremely time consuming, and a small analysis such as this may re- 
quire dozens or even scores of hours to complete. 

When the true values for the communalities were placed in the 
diagonals of the correlation matrix, no difference was found between 
the results obtained by the two methods. After rotation the factor 
loadings found by the two techniques were almost identical with the 
original factor loadings. Again it was found, however, that the Ho- 
telling method was much more time consuming. Thus it would seem 
that could a method be devised for determining the communalities, 
either method could be used to extract factors which after rotation 
would be entirely accurate. 


Summary 


Using an hypothetical three factor matrix and its corresponding 
correlation matrix, Thurstone’s method of multiple factors and Ho- 
telling’s method of principal components were applied, and the re- 
sults compared with the postulated factorial matrix. It was found: 
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1. That the third factor residuals, which ideally should be zero, 
averaged .023 by the Thurstone technique and .096 by the Hotelling 
technique. 

2. That neither set of unrotated factor loadings corresponded to 
the postulated factor loadings. 

3. That the average difference between the true loadings and 
the Thurstone factor loadings after rotation was .025, while for the 
Hotelling loadings this difference was .054. 

4. That the average difference between the true communalities 
and the communalities found by the Thurstone method was .056, 
while this difference for the Hotelling communalities was .125. 

5. That when the true communalities were used in the diagonal 
cells, the two methods were equally accurate, as was to be expected, 
but that the Hotelling method required much more time. 


TABLE I 


Hypothetical Factorial Matrix 


I II III he 

1. .900 81 
2. .800 .64 
3. -700 49 
4. .800 64 
5. .900 81 
6. .600 36 
t 424 424 36 
8. .566 .566 .64 
9. 495 495 49 
10. .520 .520 .520 81 

TABLE II 


Hypothetical Correlation Matrix 
1 2 3 4 5 6 7 8 9 


-720 
000 .000 .560 


000 .000 .000 .000 .540 

000 .000 .297 .3839 .382 .255 

509 .4538 .3896 .453 .000 .000 .240 

445 .896 .000 .000 .445 .297 .210 .280 

468 .416 .864 .416 .468 .312 .441 .588 .514 


SP Se SS 


i 





SOMO AR oP wry 


_ 


SPP Sere one 


I al ad A eh dy de oe 


— 





I 
-539 
510 
410 
439 
447 
366 
491 
661 
584 
862 


I 
.860 
.819 
014 
028 
010 


—.008 
—.011 


587 
516 
489 


I 
.614 
576 
421 
457 
448 
347 
512 
731 
648 
895 


TABLE III 


II 
—.670 
—.640 

310 

321 

345 

304 

407 
—.223 
—.193 

064 


TABLE IV 


TABLE V 


II 


-666 
658 
—.334 
—.345 
—.485 
—.A39 
—.534 
221 
123 
—.097 
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III 


-035 
023 
000 
004 
808 
-645 
466 
—.014 
512 
537 





Unrotated Factor Loadings (Thurstone) 


III 
—.048 
—.034 

514 
539 
—.571 
—.436 
—.022 
397 
—.389 
—.046 


h 


-742 
671 
528 
586 
645 
416 
407 
644 
5380 
-749 


Rotated Factor Loadings (Thurstone) 


II 
—.033 
—.025 

-726 

-766 
—.004 

029 

434 

548 
—.038 

A771 


h2 


742 
672 
527 
587 
653 
417 
406 
645 
530 
-749 


Unrotated Factor Loadings (Hotelling) 


Ill 
—.083 
—.082 

626 

.642 
—.592 
—.543 

035 

396 
—.461 
—.015 


h2 


827 
771 
681 
-740 
-786 
608 
549 
-740 
648 
811 








PPS Oe SS 


i 
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TABLE VI 
Rotated Factor Loadings (Hotelling) 
I II III h2 
1. 906 —.055 .063 .828 
2. .874 —.053 .046 -769 
8. .034 824 —.021 681 
4, .050 .859 —.006 -740 
5. —.060 .000 .885 -787 
6. —.094 —.035 Ry af: .608 
a —.054 519 525 548 
8. -653 .558 .029 -739 
9. 527 —.085 -605 .651 
10. .527 ATT 552 .810 
TABLE VII 


Comparison of Computed Thurstone and Hotelling Communalities with 
True Communalities 


True Thurstone Hotelling 
Communailities Communalities Difference Communalities Difference 

.810 -742 .068 .828 .018 
.640 671 .031 -769 129 
490 528 .038 .681 .281 
.640 586 .054 -740 .100 
.810 .645 .165 -787 .023 
.360 .416 .056 .608 .248 
.360 407 .047 548 188 
-640 .644 .004 789 .099 
490 .530 -040 .651 161 
.810 -749 .061 .810 .000 


Ave. =.056 Ave. = .125 
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EXCITATION OF A CHAIN OF NEURONES 


A. S. HOUSEHOLDER 
The University of Chicago 


When the first of a uniform chain of neurones is excited with 
a stimulus intensity S,, it is found that, according to the nature of 


the chain, either the excitation is always transmitted with decreasing 
intensity, or else the intensity for succeeding neurones of the chain 
approaches a fixed finite value, characteristic of the chain and inde- 
pendent of S,, whenever S, exceeds a certain “threshold” value. 


Consider an afferent neurone acted upon for an interval of time 
by a stimulus of constant intensity S, and developing a substance or 
“state” « (or « — 7) which serves to excite a second neurone. The 
amount of « (or e — 7) is regarded as a measure of the stimulus-in- 
tensity thus acting upon this second neurone, and it may therefore 
be denoted by S. For a particular neurone S will be a particular func- 
tion of S: 


(1) S=q(S). 


Rashevsky has considered various forms for this function [Mathe- 
matical Biophysics, Part III], one form being 


(2) S=P= (1 — e298), 


The significance of the constants need not concern us here, except for 
h, which is the threshold of the neurone, viz., that value which S must 
surpass if excitation is to take place at all. It is understood that the 
constants are all positive. 

Suppose that after suitable transformations are made S and S 
have the same physical dimensionality and are measured in the same 
units. Suppose, in fact, that this unit is h: 


(3) h=1. 


Then, in Rashevsky’s discussion, the function ~(S) always has the 
following properties: 


(i) g(1) = 90; 
(ii) yp (S) is properly monotonically increasing for S = 1; 
— 
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(iii) gy (S) is properly monotonically decreasing for S = 1; 
(iv) a = q'(o) <1. 

In fact, in the case of (2), we have 

(iv’) a=y(w)<o, 


which implies that a’ = 0. We shall not assume (iv’), however, unless 
explicitly stated. 

Suppose, now, that we have a chain of neurones, N,, Nz, --- , 
Nn, +++ , where N, excites N., N. excites N;, etc. Suppose, further, 
that these neurones are all characterized by the same “transmission” 
function »(S), so that if S; is the intensity acting upon Nj, then the 
intensity acting upon N;j., is 


(4) Sia = (Si) (t= 1,---,,+-). 
Fixing our attention upon any neurone, we may speak of S as the 
“acting intensity,” m(S) as the “transmitted intensity,” and 

(5) A(S) =g@(S) —S 

as the “intensity increment.” But apart from terminology, the recur- 
sion formulas (4) serve to define each Sj,» as some function of S;: 
(6) Siw = gr(Si) , 


where the function q, is simply @ itself. We wish to call attention to 
a few simple properties of the set of functions q, so defined. We first 
fix our attention upon the function A(S). 

On differentiating (5) we obtain 


(7) A'(S) =¢'(S) —1. 

Hence A’(S), like y’(S), is decreasing, and the graph of A(S) is con- 
cave downward. Moreove A’(S) approaches a’ — 1 asymptotically 
from above. If g’(1) > 1, then A’(S) has a single zero, Smaz > 1: 

(8) A’ (Smaz) =0, 

and A(S) a single maximum, 

(9) Amaze = A(Smez) « 

Otherwise, define 

(8’) Smaz =1, 

(9’) Amaz = —1. 


In either case A(S) decreases monotonically from Amaz for S = Smmaz 
to —o for S = «. Moreover A is everywhere negative or zero unless 


a 
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Amaz > 1, but if Amaze > 1 then there is a single interval, say S’ < S 
< S”, over which A is positive. 
We now prove 


LEMMA 1. Given any initial intensity 8S, > 1, the increments 
A(S:), A(S2), ---, are either all zero and S, = S. = ---, or else they 
are all different from zero and have the same sign, at least so long 
as the transmitted intensity exceeds the threshold unity. 

For starting with any acting intensity S, the increment is A(S) 
and the transmitted intensity is S + 4A(S). Then the succeeding in- 
crement is 


A[S+ A(S)] = g[S+ 4(S)] — [S+ A(S)] 
= ~[S+ A(S)] —¢@(S) , 
(10) A[S+ 4(S)] = 4(S) o'[S+uA(S)] (0SuS1) 
by the Theorem of the Mean. Since 
p[S+uA(S)] >0 


by (ii), it follows that 4(S) and A[S + A(S)] are both zero or else 
both are different from zero and have like signs. An obvious induc- 
tion completes the proof. 


LEMMA 2. If Anax = 0, then for any S > Snax; 
S+ A(S) > Smaz - 
If we write this inequality in the form 
S—Smaz > | A(S) |, 


it is geometrically evident, since the graph of A(S) is here concave 
downward and falling with a slope numerically less than unity. 


THEOREM. Let the sequence S, of transmitted intensities of ex- 
citation be defined by (4) with a given initial intensity S,, the func- 
tion y(S) having the properties (i) to (iv). Let the increment A(S) 
(or the decrement —A(S)) be defined by (5). The following cases 
may be distinguished: 

A. Ama < 0. Then the sequence S,; of intensities is decreasing 
until some §, falls below the threshold, and the decrement A(S,) de- 
creases as i increases as long as Si; = Smx. Since the minimum 
decrement is | Amax | , the number of excited neurones of the chain 
is certainly not greater than 


1+ (S, — 1)/| Anas | » 
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B. Amx = 0. Then if 
1. S.:< Snax, the sequence S, of intensities is decreasing with in- 
creasing decrements until some §S, falls below the threshold; 
2. S: = Brees ’ every Si = Sax 3 
3. S: >Snax, the sequence S, of intensities is decreasing and has the 
limit Smax- 
C. Amx > 0. Then A(S) has two zeros, S’ < Smax and S” > 
Smaxe If 
1. S, < S’, then the sequence §S, of intensities is decreasing with in- 
creasing decrements until some §, falls below the threshold; 
2. S.=S’, theneveryS; =S'; 
3. S’'<S, <8", then the sequence §, of intensities is monotonically 
increasing with the limit S” ; 
4. S,=S"’ ,then every 8, = S’; 
5. S, >S", then the sequence 8; of intensities decreases monotonic- 
ally with the limit S” . 
Lemma 2 is required for part B3 of the theorem; Lemma I for 
parts B2 and C. Otherwise the proof requires no further discussion. 
Let us emphasize the fact that it is here assumed that the con- 
stant initial stimulus S, is supposed to have been acting long enough 
for all neurones of the chain to have reached their stable asymptotic 
state. It is evident that previous to this time each S; will be some 
function S;(t) of the time. These functions we shall not discuss here. 
The following corollaries are almost immediate. 


COROLLARY 1. If ~(S) has a finite asymptotic value a, then, in 
case A of the Theorem, for however large an initial stimulus S, , the 
number of excited neurones of the chain is always less than 


2+ (a—1)/ | Amax |. 
In fact, if we define sequentially 
a,.—-@ 
Gin. = p(a) (i = 1, 2,---,n) 


where a, > 1 but ans, = 1, then n is the maximum possible number of 


excited neurones. Furthermore a, is the maximum possible intensity 
transmitted by neurone N,. Finally the sequences 


i(S;), y2(S:), eich »Pa(Si) 
and 


gi (81), po’ (Si), ++: »Pu (S,) 








“Seiad 


1S 
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are both decreasing for S: > Simax , and the functions g,(S),---, @n(S) 
all possess the properties (ii) and (iii). © 


COROLLARY 2. In case B of the Theorem the infinite sequence of 
functions 


pi(S),q2(S),--- 
has as its limit the function f(S) defined by 
f(S) =0, 1=8<S..:, 
£(S) =Smax , S 2 Bia, 


and the approach is monotonic for every S. 
COROLLARY 3. In case C of the Theorem the infinite sequence of 


functions 
¢:(S), g2(S),-:- 


has as its limits the function £(S) defined by 
f(S) =0, :s8<7,, 
f(S) =S’, s2s', 
and the approach is monotonic for every S. 


The dependence of the excitation upon the number of neurones 
was considered previously, for the case when ¢ is linear, by Herbert 
Landahl in a paper, “The relation between the intensity of excitation 
and the number of neurones traversed,” to appear in the September 
issue of Psychometrika. 
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EXTENSIONS OF FACTORIAL SOLUTIONS 


HARRY H. HARMAN 
The University of Chicago 


A method is developed for extending any type of factor solu- 
tion to new tests. The theoretical basis for this approximating 
scheme is thoroughly investigated, and then a simplification in the 
technique is introduced for practical purposes. An example is pre- 
sented which illustrates the procedure of extending a factor solution 
to three new tests simultaneously. 


1. Introduction 

The writer was asked by a psychologist who was interpreting a 
factor solution whether it was possible to obtain the correlation of a 
test that was not previously in the test battery, with one of the fac- 
tors. I thought the problem had practical value and sufficient statisti- 
cal interest so I undertook to solve it. While putting the final touches 
on my method of solution, my judgment as to the desirability of an 
answer to this problem was corroborated. A similar article by Paul 
S. Dwyer appeared in this journal.* The method presented in the 
present article is quite different from Dwyer’s method, although each 
has its basis in the least square theory. 

Before proceeding to an analysis of the problem, it might be well 
to note when such a problem may arise. It may happen that at the 
end of a testing program in which a factorial solution was obtained, 
a new test is given to the same subjects. Then the factor solution can 
be extended to include this test by the method to be described. All 
that need be done is to obtain the intercorrelations of this new test 
with the old tests, but no new factor analysis need be performed. An- 
other occurrence of this problem will be in the analysis of a battery 
of tests involving some general intelligence tests. General intelli- 
gence tests complicate a factorial solution in a geometric or algebraic 
sense. The complexity of such a test is equal to the number of fac- 
tors, while “in a simple structure every trait is of complexity less than 
r [number of factors].”+ Such tests of great complexity can be omit- 
ted from the factor analysis and then the factor solution can be ex- 

* Dwyer, P. S., “The Determination of the Factor Loadings of a Given Test 


from the Known Factor Loadings of Other Tests,” Psychometrika, 1987, 2, pp. 


173-178. 
+ Thurstone, L. L., The Vectors of Mind. Chicago: The University of Chi- 


cago Press, 1935, p. 155. 
aon | 
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tended to include all the tests. The descriptions of the general intelli- 
gence tests may then be in terms of every factor. 

The technique of extending a factorial solution which will be de- 
scribed is not limited to any particular type of factor analysis. As- 
sume that a factor analysis has been made of the n—1 tests 1, 2,-:-, 
n—1 and that then the correlation of test x with a factor F is desired. 
The necessary statistics are the intercorrelations of the n tests 1, 2, 

- , n—I1,n and the correlations of the original n—1 tests with F. 
The determinant of intercorrelations between tests, augmented in the 
first row and column by the correlations of the tests with the factor 
F, is written as follows: 


| 1 Tr, Tre ++ Trn-1 « 
| Tip 1 749° 99 Ti n-1 Tin 
9 Tor N21 1 pt T2,n-1 Ton 
Tn-1,P Tn-1,1 Tn-1,2°°° 1 Tr-1,n 

% Tr. Tao °°° Tana 1 








where x is the desired 7, and all other correlations are known. Whence 
A is a quadratic function in z. , 

We propose to solve this problem under the condition that the 
value of x shall make the multiple correlation coefficient Ryriis...n) a 
minimum.* In the following section we shall present the theory bas- 
ic to this assumption which yields the desired solution. Then we 
shall modify the procedure to obtain a practical solution, and finally, 
in the last section, we shall present an example of actual data illus- 
trating the method. 


2. Theoretical Treatment 


In order to simplify the following work, we make one definition. 
Let Ai;(7, 7 = F, 1, 2,--- , m) be the cofactor of 7;; in 4. Thus, Arr 
is the cofactor of the element in the first row and column of A and 
Ap» is the cofactor of & = Trn. 

The multiple correlation coefficient, in terms of determinants, is 
given by the following formulat 


* A side condition of some sort is required to solve the problem. Dwyer (op. 
cit.) assumes that the observed correlations of the new test with the original 
tests are equal to the respective correlations obtained from the factor pattern, 
without any residuals. It is thus seen that his technique differs from the one in 
this paper from original assumptions. 

+ For a discussion of the solution of partial and multiple correlations by 
means of determinants see Holzinger, K. J., Statistical Methods for Students in 
Education. Boston: Ginn and Company, 1928, pp. 312-318. 
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1 
(1) Rp (12---n) — a VArr —A. 
V Arr 

This RF is a function of x since it involves A which is a function of x. 
Now the problem is to find a value of x which will minimize R. Ap- 
plying the methods of maxima and minima of functions,* we first 
differentiate R with respect to x, set this derivative equal to zero and 
solve for x, and then put this value of x into the second derivative of 
R to see if it is positive or negative. 

Differentiating R with respect to x, we have 


dA 
dR ~ axe 


(2) tena es = 
aX 2V Arr VArr—A 








Now the derivative of A with respect to x is given by+ 


0 Tr, Tro °° Trn-1 & | 1 Yr, Vr. °° ranail | 

O 1 Me 2** Tyna Tin | Tyr 1 Tie °° Tin 9) 
dA _ |0 Too 1 ee* Tana Ton | Yor Yor Ll e+ %on10) 
ge i ee ea sg ae a ee a 

0 %n-1,1 %n-r2°°* 1 Toa |Tn-1,F Traatraaie'' 1 0 

1 Tar Tne °** Tn,n-1 1 | o Tee) Tage Tn n-1 0 | 














which, on expanding the determinants, reduces to 


dA 
(3) a = Anr + Arn = 2A rn 


since A is symmetrical. Substituting the value of = from (3) into 


(2), we have 
(4) qk = —Arm ‘ 
dt Arr VArr—A 


The derivative of R with respect to x is thus seen to vanish when 





(5) Arn — 0 ’ 
which is a linear equation in x and thus has only one root. 


* A general discussion of the maxima and minima of functions of a single 
independent variable is presented in Granville, W. A., Smith, P. F., and Longley, 
W. R., Elements of the Differential and Integral Calculus. Boston: Ginn and 
Company, 1934, pp. 182-184. 

+ The derivative of a determinant of the nth order is the sum of the n deter- 
minants obtained by differentiating each column, in turn, and leaving all the 
remaining columns unchanged. 
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We may set 
(6) Arp, = —ax-+b, 


where a is the complementary minor obtained by deleting the first and 
last rows and the first and last columns from A, and hence is positive. 
The coefficient of x is always negative as can readily be seen. If the 
order of A is odd then the sign attached to Ar, is plus, the order of 
Ap, is even, and the sign attached to the minor of x in Ap, is minus, 
so that the sign attached to a is minus. If the order of A is even then 
the sign attached to Ap, is minus, the order of Ar» is odd, and the sign 
attached to the minor of x in A,, is plus, so that again Ar, = —ax 


+ b. 
If we put the value of Ar, from (6) into (4), we finally have 


aR ___aa—b 
dz App VArr—A 





(7) 


Differentiating this expression, we obtain 


dA 
™ ER * 2a (Arr — 4) — (—ax-+ b) = 


dt 2/Anv (Arr aye 


Now we have to inspect the second derivative to see whether it is posi- 
tive or negative for x = b/a, which value makes the first derivative 
vanish. It is sufficient for this purpose to merely look at the numera- 
tor in (8). For this value of x the second term is zero. It has been 
remarked in the preceding paragraph that a is positive. Further- 
more, (Arr — A) > 0 or else the multiple correlation coefficient 
would be imaginary as is evident from formula (1).* Hence the sec- 
ond derivative is positive and R is a minimum for x = b/a. 

A multiple correlation coefficient is increased on the addition of 
a test whether this new test correlates positively or negatively with 
the preceding tests. Hence to ask for the multiple correlation co- 
efficient to be a minimum on the addition of a new test is the same as 
to ask for the multiple correlation to remain the same with the new 
test included as it was without this test. This is actually the require- 
ment we shall make in a practical problem. Instead of inquiring di- 
rectly for the minimum value of the multiple correlation coefficient 
we shall set the multiple correlation of the factor with respect to all 
the tests equal to the multiple correlation coefficient of the factor with 





*If App — 4 = 0 then Rp... ») = 0 and there is no maximum or mini- 
mum. 








be, bee 
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respect to the tests in the original analysis, i.e., with the new test ex- 
cluded. 


3. Example 


To illustrate the technique of extending a factor solution, we con- 
sider an example of 355 cases who have taken forty-one tests and 
three group intelligence tests, namely Otis (advanced), Kuhlmann- 
Anderson, and Morgan, which we shall designate by O, K, and M, 
respectively.* The intercorrelations of all tests are known and the 
factor pattern (or structure, since the factors are uncorrelated) of 
all but the intelligence tests is also known.+ Then the problem is to 
find the factor weights, or correlations with the factors, of the three 
intelligence tests. The extension of the original factor pattern to all 
three tests O, K, and M can be done simultaneously. 

Since the calculation of regression weights is extremely labori- 
ous for forty-one variables, we first reduce the number of variables 
without seriously lowering the multiple correlation coefficients. This 
reduction in the number of variables is effected by adding the tests 
that best measure each factor and dividing this sum by the standard 
deviation of the sum to get new standardized composite tests. In this 
example the groups of tests (1, 6, 20, 26, 27, 28, 29, 34, 35, 36a, 74, 75), 
(2, 3-4, 8), (11, 13, 17, 18, 22, 25b, 77, 78, 80, 81, 82, 94), (53, 
54, 55, 60), and (61, 62, 67, 68) were reduced to the composite tests 


TABLE I 
Intercorrelations of Tests 











A E V | T D 
570 
562 | .697 


593 579 611 
328 619 577 421 


550 733 891 585 .614 
591 .672 839 601 596 
507 722 824 617 655 




















BAO OCA <p 





* The data of this example are taken from Preliminary Report on Spearman- 
Holzinger Unitary Trait Study, No. 9. Prepared at the Statistical Laboratory, 
Department of Education, The University of Chicago, 1936, Tables 3 and 8. The 
factor weights of the intelligence tests as given in Table 8 were not obtained 
directly by bi-factor analysis (the procedure is described on p. 3 of Report 9). 

+ By a factor pattern is meant the set of linear equations expressing the 
tests in terms of the factors, and by a factor structure is meant the table of cor- 
relations of the tests with the facors. The elements of a factor structure are 
identical with the coefficients in the factor pattern if, and only if, the factors are 
orthogonal or uncorrelated. 
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A, E, V, T, and D, respectively.* The remaining six tests were not 
employed since they did not measure any group factor significantly. 

The intercorrelations of the five composite tests together with 
their correlations with the three intelligence tests are presented in 
Table I. The reduced factor pattern is given in Table II, where the 
general, mental speed, motor speed, spatial, verbality, attention, and 
mechanical ability factors are represented by u, a, m, e, v, t, and d, 
respectively. 


TABLE II 
Reduced Factor Pattern 











Test | u | a | m | e v t d 
A | 678 | 625 | .160 | 

E | .845 | 881 

V | 835 484 

T | .766 | All 

D | 599 | 536 














The regression weights will be computed by the Doolittle meth- 
od,+ which will be modified slightly so that these weights may be ob- 
tained first without the intelligence tests and then with them. In 
Table III the preliminary work in getting the regression weights is 
outlined. Slight discrepancies may be found in Table III due to the 
fact that all the work was carried to four decimal places and then 
rounded to three for printing purposes. The x’s appearing in the 
lower half of Table III stand for the particular correlation of the in- 
telligence test and the factor of the block in which they are located. 
For example, the very last «x in the table stands for 7ya. 

From the upper half of Table III, we can write the regression 
coefficients for estimating each of the factors in terms of A, E, V, T, 
and D. For the factor u these are as follows: 


Bup.agvr = 021 
Bur.azvp = .269 
(9) Buv.acrp = .330 
Bue.avrp = .886 
Bua.evro = .106. 


* For a description of these tests see Preliminary Report on Spearman-Hol- 
zinger Unitary Trait Study, No. 1. 

For a discussion of the Doolittle method see Holzinger, K. J., Swineford, F., 
and Harman, H. H., Student Manual of Factor Analysis (Planographed). Pre- 
pared at the Statistical Laboratory, Department of Education, The University 
of Chicago, 1937, pp. 32-36. 
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TABLE III 
Work-sheet for Computing Regression Coefficients 
Al E V T D O K M U a m e v t d 
1.) .570| .562} .593) .828/ .550| .591| .507) .678) .625) 160) 2.0.) 0.) oo} oo. 
-1.|-.570 |-.562 | —.593| -.328]/-.550 | -.591 | -.507 | -.678 |—.625/-.160| 20} | | 
1. | .697| .579| .619] .733] .672| .722] .845) ....) eR lec wie. 
~.325 |-.820|-—.338] -.187/-.3814 |-.337 |-.289 | -387 |-.3857/-.091} | | ee] eee 
675 | .876| .241| .482] .419] .8385] .483] .458/-.357|_.091| .981| _...| .....] ...... 
-1. |-.558|—.3857] —.640]-.621 | -.497 | -.642}-.679| .529) .135|-—.564| | wu]. 
1. | .611| .577} .891) .889) .824] 8385) 2...) 00] oo | SE ae 
-.316 |—.833] —.184]-.309 |-.3832 |-—.285 |-.381 |-.351/_ 999) 0) oo} | o.. 
-.210 |—.134| —.241]-.234 |-.187 |-.241]-.256| .199) 951\-212|) _...] ...] _.... 
.475| .148] .152) .848] 820] .298] .198 |-.152/_939|-.212] .484] |... 
-1. | -301| 820|-.734 |-.673 |-.627||-.418| .3821) 982| .447|-1.02| ...| ...... 
1. .421} .585| .601} .617] .766) ....) | oo] ow. 5 
~.352 |-.195]-.826 |-.351 |-.301 | -.402 |-—371/-.9095) | | ee | cee 
~.086 |—.154/-.150 |-.120|-.155]/—.164| .127) 933/186) _...] ...] _... 
—.043 | —.046}-.105 |-.096 |-.090 |-.060 | .046) 012) .064|-.146| ......|  -..... 
519} .027| 004) .034) .072) .140|-.198)|_ 051 |-.072 |-.146| .411] ...... 
-1. |-.051}-.008 |-.066 |-.189||-.270| .382) .098| .189| .281|-.793| ...... 
1.. | 6147) .806)' 0561 2009) -2c-) on) ccc | ccccantl cece, 536 
—.108}—.181 |—.194 |—.166 | —.223 |-.205/-.053} 20.0] uo} | oe. 
—.277|—.268 |—.215 |-.277||-.294 | .228) .059/-.244] 0...) 0.0] o.. 
~.048)/-.111 |-.102 |-.095 || -.063 | .049) .012| .068 |-.155 
~.001/-.000 |-.002 |-.004/-.007| .010) .003/ .004| .008/-.021| _.... 
566] .054} .083| .113) .012| .082) .021 |-.172 |-.147|-.021| .536 
-1. |-.096| -.147 |-.199 | -.021 |-.145)-.037 | .304 | .260| .037 |-.947 
Extensions to the Intelligence Tests 
Oo| K| Mi uw | a m e v t d 
1. | Tou Toa Tom Toe r r r 
-.303 | wet le aM ee tt mab pet 
-.260 —.285 .221 .057 Se ae es a eee eee oe 
-.256 —.145 .112 .029 DO? | ee tee IP, 9 ee 
-.000 -.001 .002 -000 .001 .001 “Oe to? ate 
-.005 ~-.001 —.008 —.002 .016 .014 .002 -.051 
.176 x-.806 x—-.017 x—.004 x-.063 |x-.340 | x-.001 | x-.051 
4.568 .095 .024 3 1.927 .009 .290 
-1. -5.669x —5.669x | —5.669” |-5.669x |-5.669x | —5.669x | -—5.669x 
1. Tru Ta TKm YKe Tv Txt ka 
-.849 -.401 -.370 ee ne en) eran meee es 
|-.167 —.228 177 045 WUE a Peta resteas 
|~.215 -.133 .103 .026 .143 SD ara Watt 
|-.002 —.009 .013 .003 .005 .010 —O2ZT | —  .x--- 
|-.012 | -.002 | -.012 ~.003 025 | ‘022 003 | -.079 
| .254 | 2-.773 x—.089 x—~.023 x-.016 |x-.295 | a-.024 | x.079 
| 3.040 349 .089 .064 1.158 .095 310 
| -1 | —3.9381x2 | -8.98le | -8.938lx | -—3.931x |-3.931x | -8.981x | -3.93lz 
j | by = Re ee 
1. Tuu Tua Tum Te Tuy Tut Tua 
|—.257 | —.344 -.317 ee? |, URN oF mectactilthe tas 
—278 | —.294 .229 .059 See eee rs ee ee oe 
-.187} -.124 .096 .024 133 eee te Wak are 
-.010 | -.020 .028 .007 .010 .020 wi A. paiies 
-.022 -.002 —.016 —.004 .034 .029 .004 -.107 
.246 | «—.784 x+.019 |2+.005 | x-.067 |x-.254 | x-.053 | 2-107 
3.190 —.075 —.020 Bg y 1.033 .216 435 
-1. | -4.067x | -4.067x | -4.067x | -4.067x |-4.067x —4.067x | —4.067x 
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The $ weights for the factors a, m, e, v, t, and d are similarly ob- 


tained. 
The multiple correlation coefficient for estimating a factor k 


from A, E, V, T, and D is given by 
D 

(10 R*,<aevrp) = > Bus? ix (k = U, a, M, @, V, t, d, ), 
j=A 


where x4 2vrp has been written #;,,and similarly for the other letters 
E, V, T, D. The squares of the multiple correlation coefficients thus 
obtained are as follows: 

R*y aevrp) — .8925 9 

RR, aevrp) — -7161 ’ 


RP? n(aevrp) = -0471, 
(11) RR? carvrp) = .8719 ’ 
R?yazvrp) = 5732, 
R* + anvrp) = .8269 , 
Raaevro) = .5075. 


It may be remarked that if the original factor analysis had been 
carried to completion, i.e., including the appraisals of the factors, then 
the work of the upper half of Table III would practically be done,* 
and the correlations in (11) would be known. The work of extending 
the factorial solution to the intelligence tests would then begin at this 
point. 

From the lower half of Table III, we obtain the regression co- 
efficients for estimating each of the factors in terms of six tests A, E, 
V, T, D, and O, K, or M. For each of the seven factors there are six 
beta coefficients in the extension to O, K, or M. Thus to get the seven 
factor weights for each intelligence test, 126 beta weights must be 
computed. To save space, we shall carry through the complete de- 
tails only in getting the correlation of the test O with the wu factor. 

The regression coefficients for estimating the factor u from A, E, 
V, T, D, and O are as follows: 


Buo.aevrn = —4.5680 + 5.66892 , 
Bup.aevro = .4581— .54192, 

.2856 — .0204z, 
3.5376 — 3.98112 , 
1.1483 — .945727, 
Bua.evroo .2288 — .1526v, 


where x = 7 ,. Applying formula (10) with k = wu and j ranging 
over A, E, V, T, D, O, we have 


Bur.arvoo 
(12) Buv.aerpo 


Bue.avroo 


ll 


ll 


* Except for the three columns headed O, K, and M. 
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Ry arvroo) = 4.97382 — 9.13552 + 5.66892? . 
Setting R?u;sevrp0) equal to Ru <acvro), We have 
5.66892? — 9.1355” + 4.5732 = .8925, 

or 
(13) x? — 1.611547 + .6493 = 0. 
The discriminant of equation (13) is zero, as would be expected from 
the discussion in Section 2, so that the two roots of this equation co- 
incide. Solving equation (13) we obtain 


fo = 2 = BOB. 


As a check on the preceding work, we note that the value of x 
which makes Ry :4evro0) 2 Minimum must reduce fuo.azvrp to zero and 
the remaining # values in (12) to those in (9). This is true, within 
errors of rounding, as may readily be verified. 

On the other hand, we could set Buo.azvrn equal to zero and obtain 
the required correlation x. Thus the actual computation of the mul- 
tiple correlations would be eliminated unless the complete method of 
equating the multiple correlations (with and without the new test) 
and solving a quadratic equation in x were desired as a check. Now 


setting 
Buo.arvrp = —4.5680 + 5.66892 = 0 


is equivalent to setting the negative or any multiple of this 6 equal 
to zero. Hence from Table III, we may set 


xz — .806 = 0 


and directly obtain x = .806. Thus the desired correlations may be 


read off directly from Table III. 
In a similar manner we obtain the six remaining factor loadings 


for test O, and also the factor loadings for tests K and M. These new 
approximated factor weights are presented in Table IV. 
TABLE IV 
Extended Factor Pattern 














Test | ul a m e | vw t | d 
O 806 .017 .004 .064 340 001 | .051 
K 773 089 .023 .016 295 .024 | .079 
M -784 | —.019 | —.005 .067 .254 .053 | 107 








By the method outlined in this section, we can insert a test into 
any factorial solution (orthogonal or oblique factors) and obtain 
the structure weights, or correlations of the test with the factors. 
These weights are also the coefficients in the factor pattern in the 
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case of uncorrelated factors. In case the factors are correlated, it is 
also possible to obtain the linear equation expressing the new test in 
terms of the oblique factors if desired.* 

One advantage of the method of approximating the correlation 
of a test with a factor as presented in this paper is that the correla- 
tion obtained is never greater than the correlation that would have 
arisen from the complete factor analysis including this test. This 
method of obtaining the desired correlation may then be considered 
as a conservative one.} The error in the correlation, which is bound 
to appear in any approximate scheme, is at least fixed in sign or di- 
rection. 

* See Holzinger, K. J., and Harman, H. H., “Relationships between Factors 
Obtained from Certain Analyses.” The Journal of Educational Psychology, Vol. 
XXVIII, 1937, pp. 339-341. 

+ To illustrate how the technique developed in this paper differs from Dwyer’s 


technique, the present method was applied to his example (op. cit., pp. 174-175). 
His values 7, = .3708 and r,,, —=—.7826 become Tea, — 3347 and7,,, = —.6317, 


Dwyer’s values yield zero residual correlations of the new test, while the values 
obtained by the method of this paper yield small negligible residuals as normally 
appear in a factor analysis. 
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THE APPLICATION OF MULTIPLE FACTORIAL METHODS TO 
THE STUDY OF MOTOR ABILITIES* 


CLAUDE BUXTON 
The University of Iowa 


A battery of thirteen tests of motor abilities are analyzed by 
means of the Thurstone factor methods. Interpretations of the na- 
ture of the six factors extracted are offered, with reservations. The 
high degree of specificity of motor abilities results in very narrow 
group factors. 


I. Introduction. Because most tests of fine motor ability are not 
very closely correlated, and because they lack power in predicting job 
fitness, considered either as single tests or batteries, research in the 
motor field has been neither very intensive nor very profitable. There 
has been a singular scarcity of attempts to define the concepts of abil- 
ity, performance, and capacity; to determine whether or not motor 
ability must be studied apart from intellectual ability; to determine 
whether complex motor behavior can be described in terms of ele- 
ments or basic simple activities; or to determine the effects of great 
or little motor capability upon the personality make-up of the in- 
dividual. Progress in these and similar problems seems logically to 
be dependent upon a fundamental analysis of motor performance. 

Three aspects of motor performance may be pointed out: the 
stimulus situation, the organism which acts, and the response made.t 
In the case of a worker in a factory or an S in a tapping test, the 
stimulus situation may be held relatively constant, and the response, 
which may be quantitatively scaled with a fairly high degree of ac- 
curacy, assumed to vary with the nature of the acting organism. In a 
given test situation, with motivation and like factors controlled, how 
the individual responds may be considered an indicator of: “basic” 
motor capacity, degree of skill after practice, or a combination of the 
two. These three possibilities may be termed hypotheses as to the na- 
ture of individual differences in motor performance. As in most 
studies of human abilities, we shall probably find the third to be the 

* The writer is indebted to Professors C. H. McCloy and C. E. Seashore of 
the University of Iowa, R. H. Seashore of Northwestern University, and L. L. 


Thurstone of the University of Chicago for suggestions concerning the research 


and manuscript. 
+ This description of behavior (Woodworth’s) was applied to the motor field 


by R. H. Seashore in a paper before the A. A. A. S., Indianapolis, December, 1937. 
— 
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most practicable. However, the immediate problem is to discover the 
extent to which individual differences in various kinds of performance 
vary together. Insofar as they do we can then proceed with results 
from groups of tests (whether the grouping is determined by simi- 
larity of tasks, similarity of amounts of practice, or any other fac- 
tor) to determine by various means which of the above three hypo- 
theses is the best description of the origin of individual differences. 
Factor analysis methods perhaps may indicate to us the groupings of 
motor abilities. Further factor analyses, or completely different tech- 
niques, may later be used to determine if that which is “given” in an 
individual is most important in determining his motor performance, 
or if the amount of practice he has had determines it to a greater 
degree. 

Because of the low level of intercorrelation of fine motor skills, 
individual differences have in the past been termed quite specific to 
the test. There seems to be no general motor ability which governs 
performance by a given person in many motor situations; his rank in 
the group may vary greatly from test to test. Such authors as Sea- 
shore (8) and Garfiel (5) have emphasized the probably very definite 
specificity of abilities. On the other hand, the findings of Earle (3), 
Seashore and Adams (10), Cox (2), and Farmer (4) indicate that 
coefficients may be grouped around certain “poles”, such as strength, 
rhythm, manipulative ability, steadiness, or pursuit. Such groupings 
may even in some cases be found by inspection of the correlation 
tables of those who report great specificity. The inference thus may 
be that skills sampled by such allied tests have something in common 
even though it is quite true that on the whole the level of correlation 
between varying kinds of tests is low. It is to such groupings that the 
present research draws attention. Its specific purpose is to examine 
a table of intercorrelations based upon conventional fine serial motor 
tests. Within this table it will attempt to discover any clustering of 
coefficients, broad or restricted, which indicates significant communal- 
ity of factors in performance. After this method of analysis is tried 
out we may attempt to evaluate the mode of procedure suggested 
above. 


II. Procedure. 

Setting up the test battery. In the light of previous research, the 
following hypothesis seemed to be reasonable; factors which may 
show up in motor behavior will be relatively narrow group factors. 
For this reason a test battery was set up which according to a priori 
analysis sampled several different types of action, thus providing a 
possibility for varying types of factors to appear. Furthermore, both 
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simple and relatively complex behavior was sampled, so that should 
simple tests receive significant loadings, these would be relatively 
“pure” tests and could serve as reference variables in axis rotation. 
Other general criteria for selection of tests were the following: (1) 
tests producing much fatigue were eliminated, because of the unknown 
effects of this type of variation; (2) the tasks were to some degree 
chosen for interest value, so that there were as few barriers to the 
maintenance of motivation as possible; (3) reaction-speed and choice- 
reaction tests were not used in this particular battery; (4) only tests 
for which instructions could be kept relatively simple and the proce- 
dure understandable were chosen, so as to minimize the cognitive pro- 
cesses during action; (5) only serial or repetitive action was sampled. 
Justification for such criteria of selection lies in the fact that it seems 
necessary, in the interests of clarity and economy of thinking, for an 
exploratory study to use but one general type of motor behavior. Se- 
rial action was chosen here, the tests were carefully selected, and ex- 
perimental controls enforced to make the results as easy to under- 
stand as possible. 

Apparatus. Several of the tests were constructed as modifica- 
tions of those described by other writers. (For details cf. (1), Ap- 
pendix I.) Since all are fairly well known, or the general principles 
of their construction and operation are quite simple, they will be de- 
scribed but very briefly here: 


(1, 2, 8, 4) Gains: On the four tests to be named below which showed most 
learning, viz., the rotor, cube and spool packing, and thrusting steadi- 
ness, the average of 3 or 4 of the first trials is compared with the aver- 
age of 3 or 4 of the last trials, and the percentage of improvement com- 
puted (exact number of trials for determining the averages held con- 
stant in each test, depending upon the total number of trials available 
for the computations). This admittedly is a crude set of measures; it 
is included because of the leads which it may afford. 

(5) Thrusting steadiness: S attempts to thrust a stylus into successively 
smaller holes in a brass plate; he takes 10 trials at each hole, for ac- 
curacy, from a controlled distance and at a controlled speed. 

(6) Stationary steadiness: Same as thrusting steadiness except that S holds 
the stylus in each hole for 10 seconds, unsteadiness being measured 
by the number of contacts with the plate. 

(7) Triple fre ye Three 2-in. metal discs sunk slightly below the surface 
of a board, with centers at the extremes of a 6%-in. equilateral tri- 
angle are tapped in succession with a stiff metal stylus, for speed. 

(8) Cube-packing: A low plywood box is filled with sixty-four 1 3/16-in. 
cubes, in the shortest time possible. 

(9) Spool-packing: A tray which may be emptied almost instantly is filled 
= spools a given number of times in as short an interval as pos- 
sible. 

(10) 2-bar tapping: S taps alternately between vertical metal bars, 2 inches 
apart, for speed; he uses a stiff stylus which is mounted on a light 
frame strapped to the lower arm to prevent the wrist from bending 
while the whole forearm taps. 

(11) Wrist-turn tapping: An aluminum handle is turned as frequently as 
possible through an arc of about 135 degrees between two contacts, 
with but very little resistance to the turning. 
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(12) Motility rotor: 
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ure of speed of motion in a circular path. 


(18) Pursuit-rotor: This device consists of a phonograph-like turntable car- 
rying a dime-sized metal dot, which S pursues with a stylus as it re- 


volves at constant speed. 


With the exception of the two packing tests, which were timed 
with a stop watch, every performance was objectively recorded in 
The verbal instructions, lengths of 
trials, and demonstrations by E were standardized on the basis of pre- 
vious use of these tests by the writer and others (e.g., Seashore, 9). 
A well-lighted individual testing room was used throughout the ex- 
periment, but one S being in the room with E at a time. The cycle 
method of administration of tests was used, and all tests were given 
during regular school periods, when the S’s had been excused from 


Cenco impulse-counter units. 


physical education classes. 


Subjects. School boys, 76 between the ages of 12 and 16 inclu- 


s 

ee 

13 
Reliability 

Coefficient*: .99 

No. of 
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sive, were used.* Motivation is self-provided by the tests; in addition 
an unusually satisfactory spirit prevailed among the boys, so that per- 
formance was consistently of desirable type. A more homogeneous 
age range perhaps would have been desirable, but S’s were not avail- 
able in satisfactory number except in the greater age span. 


III, Results. 

The correlation matrix. In Table I are shown the zero-order cor- 
lations which form the subject matter of this study.+ The reliabil- 
ity coefficients of the tests are very high and attest the uniformity of 
performance by the S’s. It may be noted that the correlations here 
given are quite typical of the usual fine motor skills findings, the aver- 
age 7, signs disregarded, being .15. The large number of negative co- 
efficients in Table I is due partly to the fact that most of the correla- 
tions with gains are negative — the larger the gain, the poorer must 
have been initial performance and mean score; conversely, the smaller 
the gain, the higher in the group an individual may have been in the 
beginning. Exceptions to this line of reasoning will be discussed 
later. It is obvious that if the study were to cease with inspection 
of correlation coefficients, as have so many motor skills studies in the 
past, the average r of negligible size would lead to the conclusion that 
motor abilities are not at all closely related. Even careful inspection 
of the table gives but a rough idea of possible groupings among the 
more significant coefficients. Hence it is necessary to use some means 
of reducing the zero-order matrix to more understandable simplicity. 
In this case multiple factorial methods are utilized. 

Analytic methods used. The center of gravity method of ana- 
lyzing the correlation matrix (Thurstone, 11) was used. As a means 
of reducing error in filling the diagonal cells of the correlation table 
the following technique was used: for the first factor, highest values 
in each column were inserted as communalities (11, p. 233). There- 
after the computed communalities were used, except in the few cases 
in late factors where the computed value became negative, when it 
was stepped up to zero.t A simple formula was used to determine how 
many factors should be extracted from a correlation table.§ Axes were 
rotated by plotting each variable against every other one, taken two 
at a time.|| 

* The writer is indebted to the high schools in the University of Iowa for 
fine cooperation in securing S’s, and for the use of a private testing room. 

+ All computations were rounded to four places on the basis of the four- 
place zero-order matrix. The entries in both Table I and Table II were rounded 
to two places after the analysis was completed. 

{ This method is empirically justified by Layman (7); it is partly the work 
of Professor McCloy and Miss Eleanor Metheny. 


§ Communicated to Professor McCloy by Professor Thurstone. 
| Professor Thurstone’s technique. 
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TABLE II. Corrected factor loadings 
I II III IV Vv VI h2 


1 Thrust. gain —22 —11 —.27 .29 12 .28 31 
2 Cube gain —.17 27 01 00 18 02 13 
3 Spool gain .20 26 —.05 .02 43 12 31 
4 Rotor gain .08 46 .04 .02 —.09 .02 23 
5 Thrusting 04 .09 81 .02 .03 22 83 
6 Stationary —.05 .00 82 27 25 —.02 81 
7 Triple-tap. 43 .08 15 .05 .08 59 57 
8 Cube 89 —.03 .05 —.02 .09 04 81 
9 Spool 85 .09 11 —.08 01 .09 -76 
10 2-bar —.06 Ad .00 56. .08 05 33 
11 Wrist-turn 10 31 -00 00 31 06 21 
12 Motility .28 .25 16 .29 03 57 59 
13 Rotor .00 46 12 —.01 .06 AT 44 


The factor matrix. In Table II are shown the factor loadings 
after rotation of axes. Only this table will be discussed in interpret- 
ing the findings of the study. Since only skills were studied, and skills 
ordinarily are positively correlated, the largest possible number of 
positive loadings was sought. Variable 1, gain in thrusting steadi- 
ness, proved to be the only test which could not be brought into posi- 
tive relationships with most of the primary axes. 

The heaviest loadings in factor I occur in the two packing tests, 
with the loading in triple tapping and thrusting steadiness consider- 
ably lower. At first glance this may seem like a speed factor, since 
three of the tests with the highest loadings were performed as fast 
as possible. However, the thrusting steadiness test contains very lit- 
tle of the speed element, being done at a constant, relatively slow rate. 
Furthermore, the tests which apparently call for speed and little else, 
the 2-bar tapping and motility rotor, have much smaller loadings, in 
the case of the 2-bar tapping, nearly zero. The behavior pattern 
which the more heavily loaded tests on factor I seem to demand in 
common may be described as manipulative performance or, other- 
wise stated, dexterity at speed. 

Closer examination of all the loadings in the column may give a 
clue as to what this factor may eventually turn out to be. It is pos- 
sible to conceive of a speed limit which is neurophysiological and char- 
acterizes a given restricted musculature, and which is different from 
a speed limit in which vision is the key to coordination of the hands 
in manipulating objects. Therefore this factor might be analyzed fur- 
ther by placing it in a battery with what a priori seem to be speed 
tests dependent upon eye-hand coordination in manipulating, as well 
as speed tests which require only repetitive action in a single place, or 
action not controlled by vision or audition. The test which has some- 
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thing of a loading in factor I but at present seems to have no speed 
element in it, thrusting steadiness, may be administered at different 
metronome speeds to determine more certainly whether this is true 
or not. Obviously it is not possible to name this factor at present — 
the descriptive term, manipulative ability at speed, or simply manipu- 
lative ability, explains little if anything, but it does call attention to 
the necessity for coordination of fine muscle groups, as differentiated 
from large muscle groups, e.g., in the shoulder-girdle. Further frac- 
tion is necessary to determine whether different kinds of speed (or 
other factors) are the principal determiners of performance in the 
tests more heavily loaded with I. Motion studies made with the aid 
of accurate timing devices and high speed cameras may well be the 
means of this fractionation. 

It will be remembered that as a crude measure of the amount of 
learning on four of the tests the percentage of improvement or gain 
of a few of the last trials over a few of the first trials was computed. 
When factor II is examined there is some evidence that these compu- 
tations have brought in a factor limited principally to the correlations 
with gains. The two heaviest loadings are in rotor gain and rotor to- 
tal score. The latter loading, although not in a gain, seems logical be- 
cause of the fact that each S starts from zero or near zero on the 
rotor, and total final score is positively related to the amount of im- 
provement. These loadings point to a type of improvement related 
to spatial perception. The loadings not so easily explained are 
those on wrist-turn tapping and the motility rotor, for neither of 
these tests shows much learning. However, but two loadings in the 
column reach a very significant level, they on two aspects of the same 
test, the rotor, so this factor is still quite indeterminate. 

Factor III is limited in very significant loadings to but two tests, 
and the only other loading of note is on the gain of one of these. There 
is perhaps less question, nevertheless, concerning its nature, than for 
any other factor, for the two tests are stationary and thrusting steadi- 
ness, and the loadings are relatively high. That tests of steadiness are 
apt to be more closely inter-related than many tests of motor ability 
has been shown by Seashore and Adams (10) and by Humphreys, 
Buxton, and Taylor (6). The other loading of some size, on thrusting 
steadiness gain, is negative. This is to be expected. Reversing the 
logic used in discussing the rotor gain above, we may note that in- 
stead of starting from zero or near zero and making a total score re- 
lated positively to gain, in the thrusting steadiness test an individual 
who improves a great deal usually starts with a relatively poor score; 
conversely, the person who improves little very frequently starts near 
the top of his group. Therefore the negative relationship between 
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these measures as shown in factor III is the logical one. Although eye- 
hand coordination might seem to play a part in successful perform- 
ance on these tests, the several tests which do not require much steadi- 
ness but probably require just as much eye-hand coordination have 
low loadings in III. Therefore it can be tentatively identified as a 
steadiness factor. The importance of sheer physiological tremor in 
these results cannot be evaluated at present. 

Only one loading in factor IV is of any size, that for 2-bar tap- 
ping. Although this might be considered an indication of a possible 
speed factor, or a relaxation factor, since such an ability is conducive 
to good performance on the 2-bar tapping and other tests showing 
some loadings in IV, the identification is nearly impossible. 

The only near-plausible suggestion with regard to factor V is that 
it may be concerned with a rather specific manipulative type of edu- 
cability, since it shows up to the greatest degree in the spool packing 
gain and cube packing gain. However, it is only possible to guess at 
this factor at present. This is true also of VI, which, it may be sug- 
gested, involves coordination of gross muscle groups in producing 
spatial patterns. Such an ability may be thought of as contrasted 
with factor I, where the complex coordination of fine muscle groups 
seems necessary. In the case of VI, the associated muscle groups in 
the forearm, upper arm and shoulder-girdle must control movement 
of the lower extremity, the hand, principally in one plane. This 
type of motion is called for by the motility rotor, the pursuit rotor, 
and the triple-tapping test. 

IV. Discussion. It will be seen that the factors found in this 
study must be identified with many reservations. Indeed, apart from 
one fairly satisfactory identification, the principal function of the in- 
terpretation section has been to point out leads for future testing 
which might lead to more conclusive isolations. In support of the hy- 
pothesis on which the study was founded, viz., that factors in motor 
ability are probably rather restricted group factors (limited experi- 
mentally by similarity and dissimilarity of tests, to some degree), is 
the finding, particularly in factors I and III, that large loadings, 
where they do occur, are restricted to a very few tests, and much 
smaller or zero loadings distributed to all the remainder; in no factor 
are the loadings of uniform size and distribution. 

To the criticism that such a lack of generality points to the fu- 
tility of such a method one may say: so far as the psychological prob- 
lem of the organization of motor abilities goes, broad general factors 
could not be expected to result in this study; the state of affairs is 
not so much that the methods used have failed to find the generality— 
logically the generality of abilities does not exist. Therefore the prob- 
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lem is one of discovering what the various narrow group fe ors are, 
then trying to put them into a meaningful framework. 

What is needed to make factor methods useful in evalm ing the 
three hypotheses as to the origin of individual differences in motor 
ability which were mentioned in the introduction of this paper? First 
it may be pointed out that, along the lines shown in the discussion 
above, a more restricted group of tests should be used, so that not so 
many factors are involved. Second, for each possible factor, more 
tests should be included, so that the factor may be over-determined. 
Third, the tests used should be brought nearer in content to every-day 
life situations, for it is well known that the conventional motor tests 
do not predict in any practical way the fitness of individuals for the 
tasks of industry, and it is applications of this sort for which knowl- 
edge of motor ability will be most useful. Such changes in the proce- 
dure of the present research may in the future make possible a more 
meaningful identification of factors. 
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A QUANTITATIVE DESCRIPTION OF THE STIMULUS 


J. A. LYNCH 


The Rice Institute, Houston, Texas 


A formula for the stimulus is derived. The essential postulates 
and assumptions are based on Gestalt theory. 


My present purpose is to give a comprehensive formula for the 
stimulus which lends itself to manipulation and to varied philosophical 
interpretation. I am concerned primarily with objective conditions 
which can be set in relation to human activity in the mass rather than 
with individual behavior. The only subjective factor involved in the 
general formula is intelligence which is defined statistically and dealt 
with as a constant. 

For convenience, the following symbols are introduced: 


i=n 
>> a; = A, > 
t=1 


The left-hand member is interpreted as an array of causes which 
produces the psychological effect, A. A, is the state of A correspond- 
ing to the n members of the array. 

The characteristics of the stimulus are not dependent upon the 
actual existence of A; the concern is with its objective conditions 
rather than with its concrete presence. 

The a,’s in the array are conceived as something like Gestalts 
multiplied by the sum of their “parts”. Where, 


a, = The “Whole”, and, 


a. = The sum of the parts or elements 
involved, 

@; = [A2(@, — a2) ]i 
or, 


icn 


Ay = Sai = ¥ [a,(0, —a,)],- 


This suggestion is taken, in part, from R. H. Wheeler’s first 
organismic law. 

This is of little value until the concepts of whole and parts are 
reduced to factors capable of quantitative formulation. Since the 


—95— 
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t=n 


value of these factors depends upon their occurrence in 5 a; rather 


w=1 
than in a;, it is necessary to look for measurable variables affecting 
the entire array. 

The a;’s can be described as combinations of either physical or 
sense objects (the latter including memory content, imagery, words, 
etc.) which are variously connected with the organized Gestalts. This 
probably calls for an illustration: 

Suppose, for example, that A is interpreted as an insight into 
the principle of oxidation. It will be clear later that interpretation 
as a clarified insight is only one of several possible interpretations. 

A, = an insight into this principle corresponding to 7 particu- 
lar terms of the array. 

> a; = The array of experiences, experiments, explanations, defi- 

w=1 
nitions, chemical formulas, mathematical ratios, including the learn- 
er’s own verbal and written formulations of the principle, which are 
used to clarify the principle. 

Clearly one of the factors which can vary is the total number of 
both related and unrelated objects involved in this array. This is 
the factor of complexity. Taken independently, this tends to reduce 


the value of Sa, -* Another variable is 1; this, taken independently 


(without increasing complexity) tends to increase the quantitative 


i=n 
value of > a;- Another factor to consider is the variation of the pat- 
w=1 
tern from the dominant Gestalt (R. H. Wheeler’s radial meanings). 
The latter factor is increased by the objective conditions favoring 
monotony, and particularly by the interaction of various points of 
view as in a verbal discussion between several individuals. 

These variables are not independent of each other. They can be 
grouped around two main axes of influence which are reciprocally 
related. These are the factors of concentration or intensity and dis- 
traction or complexity. (The suggestion is probably from Bechterev’s 
treatment of the relation between concentration and distraction.) 

The factors favoring intensity are as follows: 

1. A high degree of relevancy or organization within the 4@;’s, 
to reduce distracting influences to a minimum. 

2. A close analogy between the a,’s. 

3. <A large number of repetitions of situations which suggest 
the central pattern of meaning. 

4. Variation of the content or of the symbols used to convey 
the pattern to avoid polarization around elements and inhibition 


* Statement requires qualification. 


a — er 
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through fatigue. (The suggestion here is taken from the psychology 
of hypnotism, particularly as treated by I. P. Pavlov. It is also the 
process of clarifying an idea by expressing it as clearly as possible 
in a variety of ways rather than by merely repeating the same illus- 
tration or statement. Plato’s dialectical process). 

Since concentration tends to narrow the range of consciousness, 
the effects which oppose concentration can be considered as favorable 
to complexity. These factors are reciprocally related to those which 
favor intensity. They are: 

1. A large number of elements in the a;,’s. 

2. Vagueness of pattern or weak organization, allowing irrele- 
vancies which tend to generalize or spread the effects. 

3. A wide range of variation from the core pattern which also 
tends to radiate the influence over wide areas. 

These last-named factors favor a depolarization of the stimulus, 
meaning less concentration but also less inhibition; they tend to pro- 
duce a weaker but a more widely-diversified activity. 

i=n 
The stimulus, } a;, as a quantity, combines within itself both of 
=1 
these factors: Intensity which implies the tendency to call out a vig- 
orous response; and inclusiveness or breadth of content. This can be 
variously interpreted according to the value put upon A by the phi- 
losophy of education which is either implied or expressed. Only three 
of these possible interpretations are to be considered here: 

1. A =a generalized learning carried to the point of mastery. 

2. A = The total volume of a process conceived as the building 
up of associations or conversely as the breaking down of synaptic 
resistance. (The learning-process point of view). 

3. A = The total volume of activity involved (The dynamic 
point of view). 

A formula for the third criterion for the evaluation of A can be 
expressed in terms of the following factors: 

N = Complexity of the stimulus, keeping in mind the factors 
named as favoring complexity ; 

I = Intensity or polarity of the stimulus, keeping in mind also 
the factors which favor this; 

R = Resistance to learning, as synaptic resistance although not 
necessarily implying some particular nerve-path theory; and, 

M = Mental Activity (or simply activity), consciousness above 
a certain level of complexity only when activity in this sense is 
assumed to be the highest correlate of consciousness, an assumption 
which is not necessary at present. 
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The following relationships are now apparent from the stand- 
point of the activity hypothesis: 
k-NI 


oe eee (1) 


(NI) = Sa; (as a total sum). 

This formula incorporates the following distinguishable assump- 
tions: 

(a) other factors remaining constant, M varies directly with 
the intensity or polarity of the stimulus; 

(b) other factors remaining constant, M varies directly with 
the complexity or inclusiveness of the stimulus. 

(c) other factors remaining constant, M varies inversely with 
resistance to the passage of nervous impulses. 

The factors in the environment which produce (1) complexity 
and (2) polarity in the stimulus must be related closely to these con- 
cepts. All are to be considered as statistical values based primarily 
upon the conception of intelligence as a sort of a typical constant. 
This constant defines the axis of maximum effectiveness of the stimu- 
lus, determining how complex it should be on the average. This is 
not a subject of conjecture; much of the environment is already pre- 
pared as a product of cultural evolution, i.e., the accumulated, organiz- 
ing effects of the race. The problem is also subject to experimenta- 
tion; that is, various sequences of the a;’s can be tried out and 
evaluated in terms of the effect desired. This conception of the stimu- 
lus applies to the race rather than to the individual; therefore, the 
factors involved are all statistical probabilities. 

In order to get an analogous formula for the learning process 
(the breaking down of resistance, etc.), the following steps are nec- 
essary: 

Let, 

R= Rk, -&, 


where R, = the factors involved in resistance which cannot change 
with training and experience; this represents something like the re- 
ciprocal of the statistical intelligence quotient ;* and R, = The factors 
in resistance which decrease with maturity or with learning in rela- 
tion to any particular stimulus. 

R, is a constant, and R, is a variable. With respect to a particu- 
lar stimulus, the latter, R., can be reduced to zero by the learning 
process. 


* Notice the factor K as it is interpreted by diagram 2. 
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Then, 
me. 
R, = Tr’ 

in which C is a constant, and L is the learning process as defined. 

If this reasoning is sound, a substitution can now be made in 
equation J is follows: 

__k- (NI) 
MR (C/D) ’ 


Simplifying and combining constants, equation (2) results: 


K-M 
(NI) © 


But equation (2) has the disadvantage of involving two variables 
which are incapable of objective description, L and M. In order to 
get a value for L in terms of (NI) alone, it is necessary to establish 
a relationship between M and (NI) which will make it possible to 
substitute M out of equation (2). 

The general characteristics of the equation, 


f(NI) =M, 


can be derived from arguments based upon the assumptions which 
have been incorporated into the conception of the stimulus: 


L= (2) 


1. Possible Zero Points. 

By representing M on the ordinate, and (NI) on the abscissa, 
consider the possible zero points for M, i.e., the conditons which re- 
duce it to zero value: 

P, = Over-concentration (Following essentially Pavlov’s theory 
of sleep and hypnosis, monotony, etc.; this can take two forms; (1) 
repetition of a simple pattern plus variation of the symbols, as in 
hypnosis; and, (2) repetition of a simple pattern plus repetition of 
the symbols which convey it, producing radial distractions and sleep).- 

P, = Over-distraction (The condition in which each pattern is 
obscured in its turn by the pressure from other patterns. This may 
be due either (1) to lack of meaning, or to (2) over-complexity which 
obscures the core Gestalt). 


2. Possible Maximia-Minima Values. 

P, = Concentration times distraction or intensity times com- 
plexity. The inhibiting tendency of the polarizing factors is reduced 
by distraction, and the intensity of the activity is also lowered by 
distraction ; therefore, when the product of these two factors assumes 
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the highest value at which the average individual can retain the unity, 
average activity should be at a maximum. 

Since there are no other characterizing points suggested by psy- 
chological or other considerations, the curve can now be plotted.* The 
simplest relationships which will take account of the two zero points 
and the maximum value for M is the following (where a = P,, and 
b = P.): 

M = (a— (NI)][(NI) — b] = —(NI)? + (a+ 6) - (NI) —ab. 
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The above scheme gives graphic representation to the assump- 
tions, but it does not prove anything not already contained in them. 

The next step is to eliminate M in equation (2). This can be 
done by substituting (3) for M in (2), as follows: 


* The negative value of M indicated as cd can represent either over-fatigue 
or deepened hypnotic sleep from which it takes time to revive any kind of ac- 
tivity; ac can represent activity below consciousness on the over-concentration 
pole; and be can stand for lowered activity on the distraction pole. These two kinds 
of subconscious activity would evidently have different characteristics.Conscious- 
ness could also be defined arbitrarily by drawing the line at any desired level. 








lo 


(] 


(i 


ST = =— 





J. A. LYNCH 101 
_ KL— (NI) + (a+b) (NI) —ab 
L= INT) ; (4) 


Kab 
NI 


The characteristics of this curve can be determined in the fol- 
lowing manner: 





L =—K(NI) + K(a+b) — (5) 


Possible Zero Points. 

To solve for zero points or real roots, 
(1) Equate to zero and clear fraction by multiplying through by 

(NI): 

—K (NI) + K(a-+ b) (NI) — Kab =0. 
(2) Substitute in the formula for quadratic roots which is, 
—b = Vb? — 4ac 
2a ‘ 


Since a = —K, b = K(a-~- b), and c = —Kab, the substitution 
gives 


x= 








—K (a+b) = VK*(a-- 6) — 4K%ab 








aT —2K 
This gives two positive real roots when,* 
atb pees 
( = > Vab. 


Test for Maxima-minima. 
Equate the first derivative of equation (5) to zero, and solve for 
(NI), as follows: 


d(L) _ ms 
a(NI) ote aes =F - 





Ka 
NI)? 
Then, 


(NI) ==+ a = + ab. 


This tests for a maximum at -+-ab and a minimum at —ab. 


Test for Asymptotes. 
These can be found in the following manner: 


* The significance of this inequality is apparent from diagram 2. 
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(1) By inspection. When (NI) approaches zero, in equation 
(4) L approaches infinity; this establishes one of the asymptotes as 
the line, (NJ) = 0, or as the L-axis (vertical axis). 

(2) When S is the slope of the second asymptote, and B is the 
point of intersection with the vertical axis, substitute for L in equa- 
tion (5), [S(NI) + B]. This gives, 


S(NI) +B =—K(NI) + K(a+ bn Re 


Simplifying, collecting, and arranging in the order of the descending 
powers of (NI), 


(S+ K) (NI)?+ B—K(a-+ b) - (NI) + Kab=0. 
Then, set the coefficients of the two highest powers of (NJ) equal to 
zero: 

S+ K = 0, or S = —K, giving the slope of the asymptote; and, 

B— K(a+ b) = 0, or B = K( a+ b), which gives the point at 
which the second asymptote cuts the vertical axis. 

This curve can now be plotted. See the arrangement in diagram 
2. 
Comparison of the Stimulus for Activity with that for Learning: 

The graph for the L-curve, as plotted from the above character- 
istics, can be superimposed upon the same co-ordinates with the graph 
for the M-curve as shown in diagram 2, letting L share the ordinate 
with M. 

It appears from the diagram that the maximum value for L pre- 
cedes the maximum for M. This would, if proved, indicate that activ- 
ity is more compatible with complexity and perhaps with confusion 
than the learning process. This, in fact, is implied in the statement 
of the premises; but it is not so obvious until the premises are put 
together in this way. 

This seems to indicate that learning takes place best somewhere 
to the “left” of mere maximum stimulation. In fact, it can be proved 
from equation (2), which expresses the original relationship between 
M and L, that the maximum for L always precedes the maximum for 
M, for all positive values of K (K cannot be negative). 

To prove this from equation (2), which is, 


_ KM 

~ (NI)’ 

first, clear the fraction in this equation by multiplying through by 
(NI). This gives: 


L 


(NI) -L = KM. 
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DIAGRAM II* 


Second, differentiate with respect to (N/): 


dL © dM 
wn eee aay 


Now, at a maximum for L, the slope of the tangent to the L-curve, 


Be paost aL. , 
which is represented by ant is equal to zero; therefore, when L is 

* K, above, incorporates the factor described as statistical intelligence (equa- 
tion II). As expected, when that is increased the L-curve spreads out to the 
right, toward greater (complexity and intensity). 
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at a maximum, the first term 





dL 
a)» d(N1) 
is zero, and 
aM iL 
d(NI) K’ 


dM 


Since this value is positive, and since d(NT) is the tangent to the 





M-curve, the latter is still rising when L reaches its maximum value. 
This proves that the maximum value for L always precedes the maxi- 
mum for M for all positive values of K, and it can be shown by in- 
specting the graph that K cannot be negative within the limitations 
set by the problem. 

It is an important fact that this relationship between activity 
and learning does not depend upon the validity of equation (3) ; but, 
on the other hand, it helps to prove that the equation is consistent 
with the assumptions upon which equation (2) is based. 

The negative phase of the L-curve, x y z, suggests something like 
Professor Dunlap’s “Beta learnings.”* According to the above for- 
mula, these would seem to be a possibility only when fatigue or some 
of its effects act as a polarizing factor in the stimulus. In this case, 
both positive and negative learnings could result from the same stim- 
ulus as indicated by the graph. The negative aspect might take the 
form of a mild aversion. When fatigue follows practice, it might be 
expected to have this effect, even though a positive effect precedes it 
in the time sequence.+ Professor Dunlap claims that the effect of 
practice depends upon what the learner thinks and feels at the time. 
Quite so. From our point of view, this is essentially equivalent to 
saying that it depends upon the nature of the stimulus, in so far as 
thought and feeling elements enter into the complex of the stimulus.* 
The factor of fatigue could not, apparently, be incorporated in any 
positive thoughts or desires. Therefore, it would always be negative 
in its effect. A similar negative phase of the M-curve might be inter- 
preted as reflecting a similar effect upon activity. The segment, cd, 
of the M-curve would then represent the lag of M behind (NI), some- 
thing which has to be overcome by stimulation before activity can 
begin. 


ee Knight, Habits Their Making and Unmaking. Liveright, 1932. 
Pp. 16-50. 

+ Gestalts are spatio-temporal. Therefore, this statement does not necessarily 
rest upon the law of effect. R. H. Wheeler’s third and sixth laws could be inter- 
preted this way. 

* The objection to Dunlap’s theory is that he puts into the “response” ele- 
ments which functionally belong to the “stimulus.” 
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From the point of view of mastery, the interpretation of the for- 
mulas is not quite so obvious; but even here some comparisons of the 
conditions for mastery with those for learning and for activity, re- 
spectively, are more or less apparent. 

Mastery can be represented according to the first scheme by 


Sa, = A. ’ 
i=] 

which does not necessarily imply an infinite number of practice situa- 
tions, but which can be interpreted as meaning that A, = A., in val- 
ue, or very nearly so; and that A generalizes to all possible terms of 
the stimulus. This conception of mastery implies emphasis upon clar- 
ity and repeated practices; it suggests as a maximum some point ly- 
ing between V and what has been described as over-concentration, 
since this is the direction of over-learning, rather than anything on 
the confusion side of V. Therefore, it is suggested, but not proved 
conclusively as far as the relative position of W is concerned, that 
with respect to the variable quantity (NJ), W < V < U and that any 
decrease in the total quantity, (NJ), is at the expense N instead of 
I.* Since mastery is a state or condition rather than a process, a 
curve cannot be constructed for it. The process leading up to mastery 
is the learning process; this suggests that the point corresponding to 
maximum efficiency in acquiring mastered learning products is some- 
where on the left-hand slope of the L-curve. Learning as a pure pro- 
cess might be somewhat more rapid than the same process leading 
to mastery; on the other hand, it might also be attended by more 
rapid forgetting. In the latter case, some of the energy is spent in 
retaining. 


i=n i=n 
*N, =f (a,);, and 1, = f S (a, —a,) ;, following the formula introduced 
4=1 i=1 


on the first page of this paper. 
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A CONTRIBUTION TO THE MATHEMATICAL BIOPHYSICS 
OF PSYCHOPHYSICAL DISCRIMINATION 


H. D. LANDAHL 
The University of Chicago 


The properties of a simple neural mechanism are developed 
theoretically and interpreted in terms of observable variables. The 
results are discussed in connection with experimental data in cases 
for which the latter are available. 


I 


The psychological processes involved in discriminating among 
several stimuli are complex but by a series of abstractions several 
interesting phenomena may be predicted or interpreted. Let us set 
up a quantitative abstraction of a complex mechanism and systemati- 
cally investigate its properties. 

Consider the case in which two stimuli of the same modality are 
simultaneously presented for a given time, and let the magnitudes of 
a specified attribute of the physical stimuli be S, and S, such that S, 
is greater than or equal to S,. Let us assume that the excitation 
waves from the two stimuli are distinct and that they travel to a 
higher center along equivalent paths (I and II, figure 1). Then at 
the end of these paths the excitations F', and EF, are the same function 
of the magnitudes of the stimuli S, and S.. Empirically it is known 
that this function increases linearly for small values above a given 
threshold, then increases at a slower rate (1).* Let us take 


E, = K log 2! | 
(1) 
S; 
E, = K log 5- | 


as the approximate form of the relation between the excitation E at 
the higher center and the peripheral stimulus magnitude S. The 
value S must exceed the threshold h, for excitation. The general re- 


* Cf. discussion given by N. Rashevsky; and for experimental data, the 
works of Adrian and Hecht. 
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sults which follow are not greatly changed if other functions* instead 
of equation (1) are used, provided the function increases monotonical- 
ly with a decreasing slope. 

Let the paths I and II (figure 1) lead to the double synapses s, 
and s, respectively. We shall treat a double synapse as though it were 
a simple synapse. From s, let there be an excitatory fiber which leads 
to a synapse s; and an inhibitory fiber which leads to synapse s, . Sim- 
ilarly, let an excitatory and an inhibitory fiber lead from s, to s, and 
8, respectively. Let the excitatory fibers III and IV lead from s, and 
8, respectively to higher centers, so that if the excitation passes 
through III, the response FR, is made. 





FIGURE 1 


If there is no bias for either stimulus we must assume that the 
center is symmetrical. Actually there is usually a bias, as left versus 
right, but its effect can be experimentally controlled. Let us also take 
the simplest case in which the inhibitory nerve fibers have the same 
parameters as the excitatory, that is, A= B,a = b, and a and / are 


* Various relations are discussed by N. Rashevsky in chap. XXII of Mathe- 
matical Biophysics. 
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the same for both types of fibers.* When these conditions hold, then 
é; = js and e, = 7, (figure 1). The net excitations at s; and s, are 
respectively «; — 7, and «, — j;, but 


&— Js = — (3 — Ju). (2) 


Therefore only at one of the synapses can the net excitation be posi- 
tive. Hence the simultaneous presentation of the two stimuli S, and 
S, can produce only one response, FR, or R.. 

The synapses s, and s, have a threshold h,, but as long as the 
stimuli S, and S, are large enough to elicit a response, we may disre- 
gard the thresholds since only differences in the excitation factors, e, , 
and «,, will enter into later expressions. Also, the synapses s; and 
8, have a threshold h’, that is, «; — j7, must equal or exceed h’ in or- 
der for the response R, to be made. But «; — 7, will equal or exceed 
h’ provided «, — & equals or exceeds some value h. Hereafter we 
shall refer to h as if it were the threshold of synapses s; and s,. 

Using the equations developed by Rashevsky (1) we may write 
for the variation of the excitation factors at s, and s. respectively 





So = AEs, 

and (3) 
Se _ 48 — ke. 
dt 





We consider here, as throughout, that the fibers are of the simple 
type, either purely excitatory or purely inhibitory. Equations (3) 
state that the excitation factor increases with time at a rate propor- 
tional to the excitation intensity and decays at a rate proportional to 
its own concentration. The solutions of equations (3) for constant E 
are given by 

(1 — e“*) 


€ 


_ AE, 
pee g 


(4) 
(1 — e*’) 


E2 


_ AE, 
ee 


where we take «, = ¢. = 0 for t = 0 as the initial conditions. ¢ is 


* Actually it is sufficient to assume that a = b, or A/a = B/b for the fibers 
between synapses s,, 8, and s,, s,. The results in the first case are then modi- 


fied by a change in the constant factor. For notation, cf. chap. XXII, Mathe- 
matical Biophysics. 
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then the time required at this synapse only. From equations (1) and 
(4) we have 


ne So) log = (5) 


Since S, > S., the net excitation at synapse s; is positive or zero and 
proportional to (e, — h,) — (e. — hh) = e, — &, while at s, it is 
equally negative or zero and proportional to e, — e, = —(e, — &). 
We shall now introduce a concept of random distribution of dis- 
tracting stimuli at a synapse. One would expect that, at any center, 
the influence of disturbing factors, of external and internal origin, 
would not be entirely excluded. Then there would result at the center 
a random distribution of exciting and inhibiting stimuli.* Consider 
the effect of such a random influence on the synapses s, and s,. From 
the previous discussion [equation (2)] it follows that an excitation at 
8, is equivalent to an inhibition at s, and vice versa. We may there- 
fore consider the randomness as at synapse s, only. If, then, the dis- 
tribution of random stimuli were Gaussian, we would have 
.. ee 
pd(e—j)’ =———e ® d(e—}j)’ (6) 
V 210 





as the probability that the random stimulus (e« — 7)’ has a value be- 
tween (« — j)’ and («e —j)'’-+ d(e—7)’. Instead of the above equa- 
tion, let us take an approximation to it which can be handled analyti- 
cally. We take the equation 


pd(e— jy =Serewd(e—j)’,  k>0 (7) 
where 


[ipate—iy'=2 [rae =1 (8) 


and k = 1.14/c if o is the standard deviation in equation (6). 

It should be emphasized that the use of the absolute value rather 
than the square of (e — 7)’ in the exponent of the distribution func- 
tion leads to a curve which does not drop off rapidly enough for large 
values of (e — 7)’. Asa result the integral curve does not approach 
its limits rapidly enough. If such a discrepancy occurs in experiment- 
al data, it could be interpreted as evidence that the use of the normal 
distribution function would give results closely corresponding to fact. 

Positive values of (e — 7)’ act to produce an effective increase 


* The random fluctuation of the threshold is partly due to the variation in 
the metabolic activity within the nerve cell. For simpler cases of discrimination 
it is likely that this part of the fluctuation is dominant. 
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in the initial value of S, — S., strictly speaking «, — e, ; whereas a 
sufficiently large negative value of (« — 7)’ would cause S, to be re- 
sponded to as if actually less than S,. Assuming, first, that the 
threshold h is negligible, the wrong response R, will be produced 
whenever (e — 7)’ is less than —(e, — e.). The assumption that h 
is negligible is equivalent to assuming that the synapses s, and s, are 
excited equally by an amount h’ (equivalent to h) from some con- 
trolling center. But the probability, P,, that (e — 7)’ be less than 
a particular value —(«, — e.) is given by integrating equation (7) 
from negative infinity to —(e, — &.), or what is equivalent because 
of symmetry, 


Py= | pd(e—j)'’ = Werere) <i (9) 
(€3-€2) 
or 
k(e,—€) +log2P,=0. (10) 


But it follows from the above discussion that P, is the probability 
that a response is made to S, as if S. were greater than S,, whereas 
S, > Sz; or we may say that P,, is the probability of response R, to 
stimulus S,. Thus P,, is the probability of a wrong response or wrong 
choice. The definition holds whether the subject is required to choose 
the larger or the smaller stimulus. For convenience we shall consider 
the response to the greater as correct. 

Since the total excitation at synapse s, is «, — e. + (e —j)’, for 
a wrong choice to be made the net excitation at s, must be less than 
or equal to minus h, that is 


&:— & + (e—j)’S=—hw0 


or 
(e—j)' < — |e: — &| (11) 
since ¢, — €, > 0. Then using equations (5) and (10) we have 
ARK ¢y __ ¢sty log 54+ log 2P, =0. (12) 
a S2 
By a similar procedure we obtain the equation 
SSF ae log = + log 2(1 — Pc) = () (13) 


where P, is the probability of the correct response RF, . 
From equations (12) and (13) we have 
Pw + P,=1. (14) 
Let us define 
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P=2P,—1=1—2P,. (15) 


P is then the probability of a correct choice corrected for chance, that 
is, for P, = P, = 1/2, P = 0, and for P, = 1, Py = 0, P = 1. Let 
us also define 
1 

ma 
log A 
and for convenience we shall refer to D as difficulty of choice, since 
D is infinite when the stimuli are equal and D is zero* when the larg- 
er stimulus is infinitely larger than the smaller. Substituting in equa- 
tion (12) or (13) the values of P and D from equations (15) and 
(16) we have 





D= (16) 


H(i—e*) +aDlog(1—P) =0 (17) 
where 
Z-AKk. (18) 


Equation (17) defines a surface © (P,D,t) having two parameters 
H and a. This surface is somewhat similar to that discussed by L. L. 
Thurstone (2), and the variables have similar meanings but each is 
different in its definition. 


II 


Case 1. Constant t. 


For a constant time t = t, , equation (17) may be written 
Ho 
Peiait¢g F (19) 
where 


H, == (1—e*-), (20) 
Curves of equation (19) for equal increments in ¢, are represented 
in figure 2. They are initially parallel to the D-axis. It should be em- 
phasized that the D scale is quite arbitrarily defined as the reciprocal 


of the logarithm of the ratio of the stimulus magnitudes. When ¢ is 
zero, P is zero. 


Substituting equation (16) into equation (19) we have 


Hy log 53 = — log(1— P). (21) 
2 


* The minimum value of D may be very small but must be greater than zero 
since S, > h,, and S, is finite. However we shall discuss the results as if 


h,==0. 
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FIGURE 2 


If we now assign a value to P, 1 > P > 0, then 


ng = log(1-+ Si ) = constant . 
2 


2 





For small differences A S between S, and S, so that S, = S, = S, we 
may use the first term of the expansion of the logarithm and obtain 


- = constant (22) 


which is Weber’s law; and, in particular, for P. = .75 and where the 
time for judgment is long, we have for the 7. n. d., 


toe low, 2 (23) 


SS ..% 
S dH 


Case 2. Constant P. 

For a constant probability P = P,, equation (17) may be writ- 
ten 

D= B(1—e*) (24) 
where 
—H 

~ alog(1— P,) os 
Curves of equation (24) for equal increments in P, are shown in fig- 


ure 3. When D is infinite, P, = 0, and when D is zero, P = 1 unless 
t=0. 


B 0. 
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Case 3. Constant D. 


D 


FIGURE 3 





For a constant difficulty D = D,, equation (17) may be written 








ee a 
=u Se = 


(25) 
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Curves of equation (25) for equal increments in D, are shown in 
figure 4. 

The value P is defined as the probability of a correct response, 
corrected for chance, to a stimulus pair of difficulty D when the re- 
sponse is made after a time ¢. Experimentally it may be more con- 
venient to permit the subject to decide when to respond. In this case 
the value ¢ might be interpreted experimentally as the time ¢ or less. 
However, it may be that large values of t correspond more often, in 
this case, to wandering of attention or uncertainty, and hence the 
P(t) curve may decrease for large values of t. 


III 


Heretofore we have explicitly taken the threshold of synapses s; 
and s, as approximately zero. Let us now treat the case in which the 
threshold value is taken into account. Using the same notation as 
before, (e, — &2.) = — k“ log 2 P,, is the excitation at synapse s, due 
to the stimulus difference [equations (10), (11)], and (e — 7)’ is the 
value of the random excitation. Then the following relations can be 
verified. 

I. If (e — 7)’ < — (e, — & + h), the wrong response is 
made. 
Il. If — (e, —e,—h) > (e —j)' > — (a. — & + A), nei- 
ther S, nor S, is responded to, and we consider this to be an 
“equality response.” 
Ill. If (e — 7)’ > — (e, — & — h), the correct response is 
made. 
Now P,, is the probability of a wrong response when h ~ 0, or when 
only a correct or wrong choice is allowed. Let P,,’ and P,’ be respec- 
tively the probabilities for wrong and correct responses when the 
equality response is allowed, or when h > 0 ; and let P;’ be the prob- 
ability that an equality response is made. We then have the follow- 
ing relations which can be verified: 
Py = Pyre™<h, (26) 
and 
1=P,/+Pi+P. for (&—a) 2h. 
Case A. «, — & = 0,8, = S2, equal stimuli: 
Pn = Pa = Yee (27) 
Py -1e" (28) 


or 
kh = —log(1— P-). (29) 
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Case B. &—& <h, P< 4, P > Pi >? . 














1 
, = aati kh 
Pu=1—e*(Pe+75-) (30) 
p= em (31) 
c 4P. 
or 
PP’; = 1 ae wee 1 e-2kh 
ce 4P’. 
: (32) 
P,=1—P,— TP. e-2kh 
Case C. &, — & = 
ee = 2, 
P's = %o — 2 P*, (33) 
P’. = 1% 
Case D. &€&, — &2 >i, Prax FP. >. 
P’, = P,.(e* — e*) (34) 
P= 1— P,P (35) 
or 
igh 1—e™ 1—P, 


P’4= (e"—1) Py 


In each case we have P’,, + P+ P’, = 1. Cases A and C are 
special cases. However, we must treat cases B and D separately be- 
cause of the use of the absolute value of in equation (7). 


IV 


The results from the equations of cases A, B, C and D together 
with equation (19) may be illustrated by a comparison of tables I 
and II. Data from lifted weight experiments, using two and three 
categories of judgment, are given in table I (3). For the two cate- 
gories, the proportions of “greater” and “less” judgments are desig- 
nated by P, and P; while, for the three categories, the proportions of 
“greater,” “doubtful,” and “less” judgments are designated by P’,, 
P’,, and P’; respectively. Comparing the notation with those used in 
the equations above, we have P’; = P.; and P’, = P’., P, = P., 
P’, = P,,, P; = P. for stimuli greater than the standard; and 
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P’, = P’., P, = Pw, Pi = P’., P; = P, for stimuli less than the 
standard. 




















TABLE I TABLE II 
Experimental Data* Theoretical Values 
S P',,P, | P’, P’,, P; S P50" P’; F554 
grams | grams 
185 05 04 91 185 .03 .06 91 
05 (.95) .05 95 
190 12 18 -70 190 07 13 .80 
15 .85 11 .89 
195 15 25 .60 195 14 2F 59 
.20 80 24 76 
200 .30 (.42) .28 200 .29 42 .29 
52 .48 50 50 
205 55 35 10 205 .58 .28 14 
85 15 76 .24 
210 -70 18 12 210 79 14 07 
85 15 88 12 
215 85 09 .06 215 .89 07 .04 
.98 07 94 .06 


























The entries of table II are the theoretical values corresponding 
to the experimental values in table I, and were obtained from the 
equations by using the two values in parentheses to compute the two 
constants H, and kh. A plot between log (S,/S.) and —log(1—P), 
equation (21), was made. The slope of the best fitting line deter- 
mines H,, (29.5). In this case the line was passed through the point 
in parentheses in column P,;. The remaining values of P; and P,, 
and incidentally P,,, are then determined by equation (19) together 
with equation (15) and the transformations given in the preceding 
paragraph. 

Substituting the value in parentheses in column P’; into equation 
(29) determines the value of kh, (.545). In this set of data, all the 
values of P’. are greater than one-half except one value for which the 
proportions are determined by equation (27). The other values of 
P’,, P’; and P’, are then all included in Case D and determined by 
equations (26), (34) and (385) together with the transformations 
mentioned above, and the values of P,, computed from equation (19). 

In figure 5 the corresponding columns from tables I and II are 


* According to information received from Dr. Guilford, the data were ob- 
tained for illustrative purposes, and are for a single individual who made one 
hundred judgments of each pair of stimuli. 








118 PSYCHOMETRIKA 





























1.00 
u ge 
i $ 
*e 
g ° 
ioe A , 
x 
p as 
K a: PB 
- AB = en n / 
i o: 4 S.-i PB 
F ‘fe; «7 
+Ak c a é 
Wwe 
1 l 1 rf | — 1 l L 
(4 JO 100 


Experiment ————~— 
FIGURE 5 


plotted against each other. The points are labeled on the graph, the 
primed proportions referring to the case of three categories of judg- 
ment. 

In this particular set of data the point of subjective equality is 
nearly the same as that in the physical scale so that the results may 
be compared directly with those derived from the relations above. It 
may be recalled that for convenience we restricted the discussion to 
the hypothetical case in which there is no temporal or spatial separa- 
tion between the stimuli to be compared. A more complete develop- 
ment should give quantitative relations which take into account the 
temporal sequence in presentation and the spatial separation of the 
two stimuli. 

The equations (32) and (36) form a relation between P’; and P’, 
or P’,. which can be directly checked by experiment. In figure (6) the 
predicted curve, using the value of kh which was determined by equa- 
tion (29), and the experimental points are given for the same set of 
data used above. For large values of P’, , the relation P’a(P’,) is neg- 
atively linear and for small values of P’, it is positively linear. The 
use of the Gaussian distribution, equation (6), would depress the 
curve near the extremities, producing a curve which would fit the 


data more closely. 
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Let us now explicitly introduce the time ¢ into equations (26) to 
(36) by substituting for P,, its value from equation (12). 
We then have the following relations: 


log 2 P’» = log 2 P,.,—kh 
= H —at 
a wee aay AE @ )—kh (37) 
Case A: «, — & = 0, S, = S., and the relations are independent 


of time. 
Case B: ..—a& <h,P.< % 


H t H —at 
log 2(1— P's) = log |e 4 gan? lien (38) 


ae se 
log 2 P’. = — (l—e t)—-kh (39) 
Case C:e.—e,=h 


og(l ~2P) =~ <2 (1 —e~t) (40) 
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P= % (41) 
CaseD: :—a >h, P. > % 
log 2P'.= <a (1 — e*') + log (e** — e-™) (42) 
log 2(1— P’.) ce Hie) 4b. (43) 
z aD 


The probability of a wrong response P’,,, where h > 0, decreases 
with increasing t. The probability of an equality response P’; de- 
creases with t. The probability of a correct response increases with 
the time ¢ given to make the response. From equation (26) we have 
the result that P’, is proportional to P,, for any difficulty or time. The 
relation between P’, and P, is linear for P’, > 4, i.e., when the 
choices are relatively easy. When the stimuli are equal the probabil- 
ities are independent of time. The relations are in a form which is 
convenient for direct experimental check. No additional parameters 
are introduced. All the relations are specified when the three para- 
meters H, a, and kh are given. 


VI 


Thus far we have confined the discussion to the case of discrim- 
ination between two stimuli. One would expect that the consideration 
of other cases would empirically either serve as additional experimen- 
tal check or lead to a restriction which in turn would serve to evalu- 
ate one or more of the basic parameters. 

A simple approach to the case of three stimuli would be to con- 
sider three excitation waves arriving at a center similar to figure 1; 
each wave, however, has its own initial and final path. From each of 
the three lower order synapses there are two inhibitory fibers and 
one excitatory fiber leading to the higher order synapses. A similar 
mechanism is discussed for other purposes elsewhere.* If the value 
of B for the inhibitory fibers is one-half the value of A,+ we have, in 
a manner analogous to that of the previous development, the follow- 
ing value for the net excitation at the synapse leading to the correct 
response 


, ; A 
&—Js— Jo « a (2 &: — &2 — es) (44) 


where S, > S, > S; or 8: > &. > é; and the subscripts 1, 2, and 3 re- 


* Cf. figure 54, chap. XXII, Mathematical Biophysics. ; 
+ Taking B = 1/2 A simplifies the results but is not necessary. Cf. previous 
footnote. 
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fer to the lower order synapses and 4, 5, and 6 refer to the higher 
order synapses. 

If we take the case in which an intermediate response is not al- 
lowed, the subject has two degrees of choice. As in the previous case, 
the restriction arises directly from the mechanism of the discriminat- 
ing center. The distribution function for the random stimuli is then 
a surface. However, rather than work out the results, let us take the 
case in which the stimuli are nearly equal so that we may assume the 
probability for a correct response is proportional to the small excess 
excitation at the synapse for which equation (44) holds. We then 
have, using equations (44), (1) and (4), 


P «x (&—Js— Je) 
or 
AKé A, 
P =—— (1— e") lo 
i ass 





(45) 


where & is a constant of proportionality, and P is the probability of 
a correct response corrected for chance. Similarly for q nearly equal 
stimuli presented simultaneously, we have 


— AKé —ai Sit 
for the probability corrected for chance. One might then define the 
difficulty of choice, as in equation (16), by 


D,= : ' (47) 


Si 
ee 
V 293 °° Oq 


Equation (46) should then give the }(P,D,t), surface within a re- 
stricted range and under the condition that the subject can be given 
the set to treat the stimuli as a single situation rather than as mul- 
tiple pairs. For short times the proper conditions would be satisfied. 

Again, the response to the three stimuli might be treated as the 
resultant of two separate choices. Such an approach is essential 
where fairly long time is given for the choice or where the subject is 
deliberately given, or unconsciously assumes, the set to treat the 
stimuli as multiple pairs. If S,4 > Ss > Sc, where A, B, and C refer 
to the stimuli, we may make the definition that P.,s is the probability 
that the correct response A > B is made when the stimuli A and B 
are compared. We may similarly define P..;¢ and Ps.c. We then have 
as in equation (13) 








122 


or 


and similarly 


where 


Paso = 1— ze ® 
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H, 


log 2(1 — Pass) = — 


D, 


ms 


Pas = 1— se ms 


H2 


a 


Pao = 1—5e Ds 


D, = 1/log = 
B 

D, = 1/log = 
Cc 

Sp 

D, = 1/log <- 





H, = A (1 — e*4) 
a 


ate 

.-{ ow 
= —ats 
H,= = (1 — es) 








(48) 


(49) 


(50) 


(51) 


(52) 


Let us suppose that a pair of choices are made to obtain the cor- 
rect response. We have the following possible manners and orders 


of choice. 


1. 


AWW AP PD 
VV ¥¥ VV VV 
Be QAR BW AN 


BA AA Be BB 
V¥ VV VY V¥ 
Ay WER BQ WA 


BA AA BW WY 


Y¥ VV VY 
QW RW WA BO 
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In the first column A and B are compared first, and the one 
judged larger is then compared with C. In the second column A and 
C are first compared, while in the third column B and C are first com- 
pared. We assume that each of the three combinations is equally 
probable. The combinations 1, 5, 10, and 12 give the correct result 
that A is the largest, although 12 contains a false judgment that 


C>B. 
If P; is the probability of a correct response arrived at through 
the 7** combination, we then have 


Py = Ps = Pass Paso 
Py = Paso Pass (53) 
Py, = (1— Pac) Parc. 


The probability of a correct response to the three stimuli is then 


P.=5 (P, + Ps + Pro + Piz) 


Ts (4 “se 1 cya 1 se \(1 baa 
= = —m = EC 1 msm C 2 — = 1 —-C 3 
al ( 2 \ 2 oe 2 ( 2 


oe. 1 4% 
<2 (1—~e ®) 
5 ( 2 
or 
a ae cee mm | A.B 1  * 
P, = Junwe DoH —@e Ds —e 1 Dd cage Ds Dean ag De Dee, 
2 3 TS Tp 12 
(54) 


For S, = 0, then D, = 0 and D; = 0, and we have 


Hy 


1 _»*: 
Pe=1—se ™ =Paos- (55) 


Thus the case of three stimuli reduces to that of two stimuli when the 
smallest is much smaller than the other two. If Ss; and S> are both 
negligible compared with S,, then P, = 1. 


If S, = Ss = Sc, then D,, D., and D, are large and we may re- 
tain only the first term of the expansions of the exponentials in equa- 


tion (54). We may then also assume that, on the average, t, = t = 


t, so that H, = H, = H; = Hy. Then equation (54) reduces to 
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= FS oe 
eet te 


(56) 











Hy i Sa 


_. 
3° 2 VS Se 


which is similar to equation (45) since P, is not corrected for chance. 


Then introducing P = S (P. — ; ) in the above equation 


2 
_8H - 
Pp =a Os to) (57) 
where 
D= —— (58) 
log —— 
VSe Se 


and t, is the time for each decision, or one-half the total time for re- 
sponse. The similarity between equations (45) and (57) is evident 
and the concept of D in each case is the same. Several of the equa- 
tions are in a form convenient for experimental verification. 

If A,S = S, — Sz and AS = S, — S_, the denominator in equa- 
tion (58) may be expanded, since the values S,, Sz and Sc are very 
nearly equal, and we have 


S, _148+4,8 - 
waa £ Ss 


For large values of ¢t, and for P = 14, we then have from equation 
(57) 





:. - we 


which may be compared with equation (23). 


VIll 


One would not expect the use of the Gaussian distribution [equa- 
tion (6)] to have greatly changed the results obtained. However, the 
assumption that the decay of the excitation factor proceeds at the 
same rate as the decay of the inhibition factor (i.e.,a = b) might be 
expected to lead to unnecessary restriction in the time factor. If 
a =~ b, we would have an additional parameter for the surface 
»(P,D,t) and the shape would be somewhat changed. 
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The development has been based on the assumption that the ran- 
dom excitation (e« — 7)’ has its value at the time of presentation of 
the stimuli and maintains it until withdrawal of the stimuli. If, how- 
ever, we take a more general case in which there is also a fluctuation, 
random with respect to time, we may be able to study the effect of 
drugs, motivation, distracting factors, etc., on the constants a, H and 
kh and on the constants of the time distribution function or on the 
frequency, if the simpler case of constant frequency is assumed. Such 
relations as that between difficulty and length of time for judgment 
should then be given. In this more general case the motivation may 
be either for shorter judgment time or for increased accuracy, where- 
as in the case discussed, the motivation could be only for greater ac- 
curacy. It would seem most plausible that increased motivation 
might, over a small range, increase the inhibition of distracting fac- 
tors thus decreasing k. Then, over this small range, for a given dif- 
ficulty and time, the probability for a correct response would increase 
with motivation. 

On the basis of a rather simplified mechanism we have derived 
relations between the observable proportions P,, P;, P',, Pa, P1, 
the difficulty in terms of the magnitudes of the physical stimuli 
D(S,, Sz), and the judgment time t. The relations, for which there 
are experimental data, show sufficiently good agreement (figure 5). 
The other relations conform to that which one would expect to be 
plausible from other considerations. 

This investigation was made possible by a grant from the Rocke- 
feller Foundation to the University of Chicago. 
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ERRATA 


PSYCHOMETRIKA, Vol. 3, No. 1, March, 1938 


P. 19, line 7. For “a; = 1--- n—1” read “a;,i = 1---n—1”. 


P. 20, line 2. For “the 7 + 1 square matrix” read “the n + 1 square 
matrix”’. 


P. 20, line 3 from bottom. For “is less than” read “is two less than”. 











ae 


